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Studies  using  an  existing  31  lumped  mass  model  of  the  UH-1H  helicopter  are 
performed  to  determine  the  sensitivity  of  responses  to  changes  in  the  load- 
deflection  r ‘presentation  for  the  engine  and  transmission  mounts,  landing 
gear  and  fuselage.  Simplified  techniques  are  used  to  predict  the  load- 
deflection  curve  for  the  crushing  of  a  segment  of  the  lover  fuselage  under 
impact  conditions.  The  predictions  include  elastic  behavior,  failure  load 
and  post- failure  behavior.  The  structural  segment  selected  for  analysis  and 
test  is  a  section  supported  on  four  edges,  representative  of  the  transmission 
pylon  support. 

Twelve  specimens  were  fabricated.  The  specimens  are  h6  inches  long  by  18 
inches  vide  by  6.125  inches  to  12.125  inches  deep.  The  6.125-inch-deep 
specimens  are  approximately  half  the  size  (except  for  thickness)  and  are 
varied  in  detail  design  (number  of  angles,  spacing  of  angles,  lightening 
holes).  Static  and  dynamic  tests  were  performed.  The  predicted  load- 
deflection  curves  are  compared  to  the  test  load-deflection  curves  and  show 
good  agreement  with  regard  to  peak  failure  load,  failure  point,  energy  ab¬ 
sorbed,  and  shape.  The  results  of  the  tests  shov  that  for  this  type  of 
typical  fuselage  structure,  static  tests  provide  load-deflection  data  which 
is  similar  to  data  that  can  be  obtained  from  dynamic  tests,  but  more  economi¬ 
cally. 

Program  KRASH  is  refined  to  facilitate  its  use  by  designers.  In  particular, 
the  input  data  is  reordered,  some  inputs  are  standardized  and  more  general 
load-deflection  curve  characteristics  are  possible.  The  capacity  of  the 
program  is  increased  to  80  lumped  masses,  100  internal  beams  and  120  load- 
deflection  tables.  The  refined  program  was  run  to  demonstrate  capability 
to  treat  a  three-dimensional  impact  velocity,  mass  penetration  into  an 
occupiable  space,  and  simplified  rotor  blade  contact.  Specimen  test  data 
is  also  used  to  refine  the  31  mass  UH-1H  model. 

The  analytical  techniques  developed  herein  are  presented  in  the  form  of  de¬ 
sign  charts,  nomographs,  curves,  tables  and  equations  and  form  the  basis  of  a 
structural  crashworthiness  design  manual.  'The  design  procedures  are  outlined 
in  a  step-by-step  process  including  examples 

Volume  II  contains  supporting  analytical  and  test  data  and  a  literature 
matrix  categorization.  A  description  of  refined  program  KRASH  is  provided 
which  shows  the  new  input -output  format,  a  listing,  and  sample  problems. 
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INTRODUCTION 


This  volume  contains  test  and  analysis  data  in  support  of  the  detailed 
discussions  presented  in  Volume  I.  The  60  technical  reports  and  publica¬ 
tions  reviewed  during  the  study  are  briefly  summarized  and  then  tabulated 
in  a  Literature  Survey  Subject  Index.  The  test  data  is  presented  in  the 
form  of  recorded  load,  deflection,  strain  and  accelerometer  readings 
versus  scan  (time).  The  refinements  of  program  KRASH  to  facilitiate  its 
usage  are  described.  Included  in  the  writeup  for  program  KRASH  is  the 
refined  input-output  format,  user's  guide,  listing,  and  two  sample  prob¬ 
lems.  Also  included  in  the  section  on  program  KRASH  are  the  energy  balance 
equations  and  the  sample  format  for  the  energy  data. 
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LITERATURE  SURVEY  AND  SYNOPSIS 


1.  Wittlin,  G.,  Gamon,  M.A.,  EXPERIMENTAL  PROGRAM  FOR  THE  DEVELOPMENT  OF 
IMPROVED  HELICOPTER  STRUCTURAL  CRASHWORTHINESS  ANALYTICAL  AND  DESIGN 
TECHNIQUES,  Lcckheed-Califomia  Company;  USAAMRDL  Technical  Report 
72-72A,-72B,  Eustis  Directorate,  U.S.  Army  Air  Mobility  Research  and 
Development  Laboratory,  Fort  Eustis,  Virginia,  May  1973>  AD  764985* 

This  report  presents  the  results  of  a  four-phase  study  to  develop  im¬ 
proved  rotary-wing  aircraft  structural  crashworthiness  analytical  and 
design  techniques.  A  digital  computer  program,  designated  KRASH,  was 
developed  and  shown  to  be  capable  of  accurately  predicting  responses 
during  a  crash  in  which  multidirectional  forces  are  present.  The  pro¬ 
gram  was  verified  using  controlled  test  data  obtained  on  the  same  pro¬ 
gram  via  a  combined  vertical- lateral  velocity  impact  using  a  fully 
instrumented  UH-1H  helicopter.  Parameter  studies  were  performed  to 
ascertain  the  effect  on  the  response  of  the  structure  and  the  occupants 
due  to  changes  in  structural  element  load  deflection  characteristics. 

A  consistent  design  approach  is  presented,  and  the  results  of  the  para¬ 
meter  study  are  used  to  illustrate  its  application  in  a  crash  analysis. 
The  study  also  included  a  detailed  accident  investigation,  a  literature 
survey  and  evaluation,  substructure  test,  and  analysis. 

2.  Turnbowi  J.W.,  Carroll,  D.F.,  Haley,  J.L.,  Jr.,  Robertson,  S.N.,  CRASH 
SURVIVAL  DESIGN  GUIDE,  Dynamic  Science;  USAAMRDL  Technical  Report 
71-22,  U.S.  Army  Air  Mobility  Research  and  Development  Laboratory, 

Fort  Eustis,  Virginia,  October  1971,  AD  733358. 

This  report  is  a  design  guide  that  has  been  assembled  to  provide  the 
engineer  with  an  understanding  of  the  basic  problems  associated  with 
the  development  of  crashworthy  U.S.  Army  aircraft.  Where  possible, 
solutions  to  specific  problems  are  indicated.  In  areas  in  which 
little  design  data  are  available,  only  the  general  philosophy  appro¬ 
priate  to  the  problem  solution  is  presented;  the  details  of  such  solu¬ 
tions  as  well  as  the  degree  of  crashworthiness  to  be  achieved  must  be 
left,  at  present,  to  the  ingenuity  of  the  designer. 

This  guide  presents,  in  a  condensed  form,  the  data,  design  techniques, 
and  criteria  that  are  presently  available  in  eight  areas:  (l)  air¬ 
craft  crash  kinematic  and  survival  envelope,  (2)  airframe  crashworthi¬ 
ness  design  criteria,  (3)  aircraft  seat  design  criteria  (crew  and 
troop/passenger),  (4)  restraint  system  design  criteria  (crew,  troop/ 
passenger,  and  cargo),  (5)  occupant  environment  design  criteria, 

(6)  aircraft  ancillary  equipment  stowage  design  criteria,  (7)  emergency 
escape  provisions,  and  (8)  postcrash  fire  design  criteria. 

It  is  intended  that  both  airframe  and  component  designers  and  manu¬ 
facturers  use  this  guide  to  extend  the  "region  of  survivability"  in 
aircraft  accidents  to  a  maximum  level. 
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3.  Reed,  William  H.,  Avery,  James  P.,  Fh.D,  PRINCIPLES  FOR  IMPROVING 
STRUCTURAL  CRASHWORTHINESS  FOR  STOL  AND  CTOL  AIRCRAFT,  Aviation  Safety 
Engineering  and  Research;  USAAVLABS  Technical  Report  66-39 >  U.S.  Army 
Aviation  Materiel  Laboratories,  Fort  Eustis,  Virginia,  June  1966, 

AD  637133. 

In  this  report  the  crash  behavior  of  aircraft  structures  is  investi¬ 
gated.  The  investigation  begins  with  the  definition  of  two  indices 
of  crashworthiness  of  basic  aircraft  structures  and  the  analysis  of 
the  influence  of  several  general  types  of  structural  modifications 
upon  these  two  indices.  This  analysis,  using  fundamental  principles 
of  mechanics,  contains  several  simplifying  assumptions,  which  are  ex¬ 
plained  as  they  axe  introduced. 

Design  concepts  to  improve  the  ability  of  the  "protective  container" 
to  maintain  living  space  for  occupants  during  a  crash  or  to  attenuate 
the  accelerations  experienced  by  occupants  during  a  crash  are  develop¬ 
ed  for  crash  conditions  which  are  either  primarily  longitudinal  in 
nature  or  primarily  vertical  in  nature.  Analytical  methods  are  then 
provided  to  show  how  and  when  to  apply  these  design  concepts  to  any 
particular  aircraft.  Principles  are  presented  which  are  considered  to 
be  suitable  for  use  during  design  of  new  aircraft  as  well  as  modifica¬ 
tions  of  existing  aircraft. 

The  results  are  presented  from  three  full-scale  crash  tests  of  small 
twin-engine  airplanes  which  were  conducted  as  a  part  of  this  investi¬ 
gation. 

Among  the  pertinent  conclusions  of  the  report  are:  (l)  improvements 
in  crashworthiness  can  be  achieved  via  minor  changes  in  structural 
design  or  modification  of  existing  structure,  (2)  vertical  and  longi¬ 
tudinal  impact  environments  offer  significantly  different  problems  in 
designing  for  improved  crashworthiness,  and  (3)  analysis  of  aircraft 
behavior  is  hampered  by  the  lack  of  adequate  knowledge  of  the  rela¬ 
tionships  which  apply  to  the  determination  of  the  reaction  force  which 
decelerates  the  aircraft  upon  contact  with  the  ground. 

4.  Greer,  D.L.,  et  al,  DESIGN  STUDY  AND  MODEL  STRUCTURES  TEST  PROGRAM  TO 
IMPROVE  FUSELAGE  CRASHWORTHINESS,  General  Dynamics,  Convair;  FAA  Re¬ 
port  DS-67-20,  Federal  Aviation  Administration,  Washington,  D.C., 
October  1967,  AD  666816. 

This  report  presents  the  results  of  a  study  to  evaluate  methods  to 
improve  crashworthiness  by  retaining  transport  cabin  integrity  during 
crash  impact  loadings. 

The  study  includes  analyses  of  the  effects  of  strengthening,  redistri¬ 
bution  of  bending  material,  and  incorporation  of  energy  dissipating 
features  on  the  ability  of  the  fuselage  to  provide  a  protective  shell 
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around  the  occupants.  Analytical  results  were  substantiated  by  a  test 
program. 

The  test  program  included  compression  tests  of  plate-stringer  panels 
and  drop  tests  of  representative  fuselage  structure.  Tests  were  made 
on  three  100- inch-diameter  cylindrical  sections  dropped  axially,  four 
segments  of  100-inch-diameter  cylinders  dropped  laterally,  and  a  struc¬ 
turally  comple'  nose  section  of*  a  jet  transport  dropped  in  a  10-degree 
nose-down  attiwde. 

The  requirement  for  a  plastically  deforming  structure  is  important  for 
both  axial  and  vertical  collapse  characteristics  of  a  fuselage.  Plas¬ 
tic  collapse  provides  the  most  efficient  energy-absorbing  capability 
and  also  reduces  the  possibility  of  excessive  tearing  or  complete  dis¬ 
integration  of  the  structure. 

Both  energy-absorbing  capacity  and  failure  mode  are  important  for  ver¬ 
tical  crushing  conditions.  The  recommended  manner  of  reinforcement 
for  the  fuselage  lower  frame  segments  strengthens  the  bottom  center- 
line  portion  and  the  floor  beam  to  frame  area.  No  significant  weight 
or  cost  penalty  is  involved  since  the  crash  requirement  reinforcement 
occurs  at  the  most  critical  areas  for  existing  design  conditions. 

5.  Saczalski ,  K.J.,  and  Park,  K.C.,  TRANSIENT  RESPONSE  OF  INELASTICALLY 
CONSTRAINED  RIGID  BODY  SYSTEMS,  to  appear  in  the  Journal  of  Engineer¬ 
ing  for  Industry.  A.S.M.E.  Transactions,  1974. 

An  energy  rate  balance  is  employed  to  develop  the  incremental  equations 
of  motion  for  a  shock- loaded,  inelastically  constrained,  rigid-body 
structural  system.  Lagrangian  multipliers  provide  the  coupling  mecha¬ 
nism  necessary  to  reduce  the  overall  system  of  equations  to  a  set  of 
modified  rigid-body  equations  which  include  the  nonlinear  geometric 
and  structural  material  effects.  Kinematic  material  hardening  and  a 
modified  yield  criteria  are  used.  Examples  illustrate  the  technique 
and  are  compared  with  experimental  results. 

6.  Massonnet,  C.E.,  and  Save,  M.A.,  PLASTIC  ANALYSIS  AND  DESIGN,  VOL.  I, 
BEAMS  AND  FRAMES;  Blusdell,  New  York,  1965. 

This  book  deals  with  the  analysis  and  design  of  beams  and  frames  made 
of  a  ductile  material  on  the  basis  of  the  ultimate  load.  With  the 
second  volume,  which  is  concerned  with  more  complicated  structures 
such  as  plates  and  shells,  it  aims  at  giving  a  broad  (if  not  exhaus¬ 
tive)  coverage  of  plastic  analysis  and  design  methods.  These  methods 
essentially  apply  to  mild  steel  structures,  but  may  also  be  used,  with 
adequate  caution,  for  reinforced  and  prestressed  concrete  structures. 

Presently,  the  first  volume  is  the  most  important  for  practical  appli¬ 
cations.  The  so-called  simple  plastic  theory  of  beams  and  frames  is 
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presented  In  Volume  II, which  is  now  being  finalized. 

The  material  has  been  considered  as  an  engineering  problem  and  not  as 
applied  mathematics .  Much  attention  and  extensive  treatment  have 
been  given  to  plastic  buckling  and  to  design  of  joints.  Hie  in¬ 
fluence  of  details  of  construction  (machining,  drilling,  punching, 
etc.)  is  also  considered. 

7.  Roark,  R.J.,  FORMULAS  FOR  STRESS  AND  STRAIN,  McGraw-Hill,  New  York, 

1965. 

This  book  provides  a  compact,  complete  summary  of  the  formulas,  facts, 
and  principles  pertaining  to  strength  of  materials.  It  is  intended 
primarily  as  a  reference  book  and  represents  an  attempt  to  meet  what 
is  believed  to  be  a  present  need  of  the  design  engineer. 

Presented  sure  certain  general  principles.  Included  are  brief  descrip¬ 
tions  of  analytical  and  experimented  methods  of  stress  analysis  and 
information  concerning  the  behavior  of  materied  under  stress.  The  be¬ 
havior  of  structural  elements  under  various  conditions  of  loading  is 
discussed,  and  extensive  tables  of  formulas  for  the  calculation  of 
stress,  strain,  and  strength  are  given.  Derivations  and  detedled 
explanations  are  omitted.  However,  examples  are  included  to  illus¬ 
trate  the  application  of  the  various  formulas  and  methods. 

8.  Ayre,  Robert  S.,  Shock  and  Vibration  Handbook,  Chap.  8,  Vol.  1, 
TRANSIENT  RESPONSE  TO  STEP  AND  PULSE  FUNCTIONS,  McGraw-Hill,  1961. 

Hoppmann,  W.K. ,  Shock  and  Vibration  Handbook,  Chap.  9,  Vol.  1,  EFFECTS 
OF  IMPACT  ON  STRUCTURES,  McGraw-Hill,  New  York,  1961. 

Chapter  8  deals  briefly  with  methods  of  analysis  for  obtaining  the 
response  spectrum  from  the  time  history,  and  includes  in  graphical 
form  certain  significant  spectra  for  various  regular  step-  and  pulse- 
type  excitations.  The  usual  concept  of  the  response  spectrum  is 
based  upon  the  single-degree-of-freedora  system,  usually  considered 
linear  and  undamped  although  useful  information  sometimes  can  be 
obtained  by  introducing  nonlinearity  or  damping.  The  single-degree- 
of-freedom  system  is  considered  to  be  subjected  to  the  shock  or 
transient  vibration,  and  its  response  is  determined. 

The  response  spectrum  is  a  graphical  presentation  of  a  reflected 
quantity  in  the  response  taken  with  reference  to  a  quantity  in  the 
excitation.  It  is  plotted  as  a  function  of  a  dimensionless  parameter 
that  includes  the  natural  period  of  the  responding  system  and  a  sig¬ 
nificant  period  of  the  excitation.  The  excitation  may  be  defined  in 
terms  of  various  physical  quantities,  and  the  response  spectrum  like¬ 
wise  may  depict  various  characteristics  of  the  response. 
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Chapter  9  discusses  a  particular  phenomenon  in  the  general  field  of 
shock  and  vibration  usually  referred  to  as  impact.  An  impact  occurs 
when  two  or  more  bodies  collide.  An  important  characteristic  of  an 
impact  is  the  generation  of  relatively  large  forces  at  points  of  con¬ 
tact  for  relatively  short  periods  of  time.  Such  forces  sometimes 
are  referred  to  as  impulse-type  forces. 

Three  general  classes  of  impact  are  considered  in  this  chapter:  (l) 
impact  between  spheres  or  other  rigid  bodies,  where  a  body  is  con¬ 
sidered  to  be  rigid  if  its  dimensions  are  large  relative  to  the  wave¬ 
lengths  of  the  elastic  stress  waves  in  the  body;  (2)  impact  of  a  rigid 
body  against  a  beam  or  plate  that  remains  substantially  elastic  during 
the  impact;  and  (3)  impact  involving  yielding  of  structures. 

9.  Timoshenko,  S.  Woinowskv-Krieger,  S.,  THEORY  OF  PLATES  AND  SHELLS, 
McGraw-Hill  Publishers,  New  York,  1959* 

Ihis  book  deals  with  the  three  regions  of  plate  and  shell  theory: 

(l)  thin  plates  with  small  deflections;  (2)  thin  plates  with  large 
deflections;  (3)  thick  plates.  The  book  considers  problems  with  mem¬ 
brane  stresses  and  the  case  with  clamped  edges.  Simplifications  are 
given  for  special  cases  of  deformation  to  the  shape  of  a  developable 
surface.  The  thick  plate  theory  presented  considers  the  problem  as  a 
three-dimensional  problem  of  elasticity. 

10.  Tulk,  F.D.,  BUCKLING  OF  CIRCULAR  CYLINDRICAL  SHELLS  UNDER  DYNAMICALLY 
APPLIED  AXIAL  LOADS,  UTIAS  report  160,  1972. 

A  theoretical  and  experimental  study  was  made  of  the  buckling  charac¬ 
teristics  of  perfect  and  imperfect  circular  cylindrical  shells  sub¬ 
jected  to  dynamic  axial  loading. 

The  tests  were  performed  on  a  specially  designed  dynamic  testing 
machine  which  was  capable  of  producing  controlled  ramp-type  loads 
at  rates  ranging  from  the  quasi-static  up  to  higher  than  200,000 
pounds  /second.  Ten  shell  specimens  were  tested,  including  two  near- 
perfect  shells,  seven  shells  with  axisymmetric  sinusoidal  imperfec¬ 
tions  of  a  variety  of  amplitudes  and  wavelengths,  and  one  shell  with 
quasi-random  axisymmetric  imperfections.  The  shells  were  produced 
from  a  photoelastic  epoxy  plastic  using  a  spin-casting  technique. 

The  imperfection  profiles  were  machined  into  the  shell  walls  using 
a  high-precision  hydraulic  tracing  apparatus.  For  three  of  the 
shells  with  sinusoidal  imperfections,  imperfection  profiles  were  cut 
on  the  inner  surface  alone;  while  for  the  remaining  four  shells  with 
sinusoidal  imperfections  and  for  the  shell  with  quasi-random  imper¬ 
fections,  a  special  manufacturing  procedure  was  adopted  which  pro¬ 
duced  shells  with  matching  inner  and  outer  profiles,  thus  providing 
effectively  constant  thickness  walls. 
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Experimental  data  included  dynamic  buckling  loads  (124  data  points),  | 
high-speed  photographs  of  the  buckling  mode  shapes,  and  observations  | 
of  the  dynamic  stability  of  shells  subjected  to  rapidly  applied  sub-  | 
critical  loads.  1 

) 

A  mathematical  model  is  developed  to  describe  the  dynamic  behavior  I 

of  perfect  and  imperfect  shells.  This  model  is  based  on  the  Donnell-  | 

von  Karman  compatibility  and  equilibrium  equations  and  has  a  wall  de-  | 

flection  function  incorporating  five  separate  modes  of  deflection.  ^ 

Close  agreement  between  theory  and  experiment  is  found  for  both  dy-  < 

namic  buckling  strength  and  buckling  mode  shapes.  I 

S 


11.  Stricklin,  J.E.,  et  al,  LARGE  DEFLECTION  ELASTIC -PLASTIC  DYNAMIC 
RESPONSE  OF  STIFFENED  SHELLS  OF  REVOLUTION,  TEES-RFT-72-25  and  SLA- 73- 
0128,  1972. 

This  report  presents  the  formulation  and  check-out  problems  for  a 
computer  code  DTNaI’IAS,  which  analyzes  the  large  deflection  eiestic- 
plastic  dynamic  response  of  stiffened  shells  of  revolution.  The  for¬ 
mulation  is  by  the  finite  element  method,  with  finite  differences 
being  used  for  the  evaluation  of  the  pseudo  forces  due  to  material 
and  geometric  nonlinearities .  Time  integration  is  by  the  Houbolt 
method.  The  stiffness  may  be  due  to  concentrated  or  distributed 
eccentric  rings  and  spring  supports  at  arbitrary  angles  around  the 
circumference  of  the  elements.  Check-out  problems  include  the  com¬ 
parison  of  solutions  from  DYNAPLAS  with  experimental  and  other  com¬ 
puter  calculations  for  rings,  conical  and  cylindrical  shells,  and  a 
curved  panel.  A  hypothetical  submarine  including  stiffeners  and 
missile  tube  is  studied  under  a  combination  of  hydrostatic  and  dyna¬ 
mically  applied  asymmetrical  pressure  loadings. 

12.  Becker,  H.,and  Gerard  G. ,  HANDBOOKS  OF  STRUCTURAL  STABILITY,  PARTS  I- 
V,  NACA  Technical  Notes  3781-3785,  1957. 

The  local  bv  :kling  of  stiffener  sections  and  the  buckling  of  plates 
with  angle  stiffeners  are  reviewed,  and  the  results  are  summarized 
in  charts  and  tables.  Numerical  values  of  buckling  coefficients  are 
presented  for  longitudinally  compressed  stiffener  sections  of  various 
shapes,  and  for  stiffened  cylinders  loaded  in  torsion.  Although  the 
data  presented  consists  primarily  of  elastic-buckling  coefficients, 
the  effects  of  plasticity  are  discussed  for  a  few  special  cases. 

13.  Skogh,  J.,  Stem,  P. ,  POSTBUCKLING  BEHAVIOR  OF  A  SECTION  REPRESENTA¬ 
TIVE  OF  THE  B-l  AFT  INTERMEDIATE  FUSELAGE,  Lockheed  Palto  Alto  Re¬ 
search  Laboratory,  AFFDL-TR-73-t>3*Air  Force  Flight  Dynamics  Labora¬ 
tory,  Air  Force  Systems  Command,  Wright-Patterson  Air  Force  Base, 

Ohio,  May  1973. 
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A  spction  of  the  B-l  aft  intermediate  fuselage  consisting  of  a  com¬ 
bination  of  flat  and  curved  panels  is  analyzed  for  postbuckling 
strength  under  a. combination  of  torque  and  axial  loading.  The  analy¬ 
sis,  which  extends  to  load  levels  about  ten  times  the  load  that  pro¬ 
duces  the  first  buckle,  was  carried  out  rigorously  by  the  use  of  the 
finite-difference  computer  code  STAGS.  The  results  show  that  the 
fuselage  section  does  not  collapse  at  the  applied  ultimate  load. 

Shear  stiffness  values  as  a  function  of  the  applied  load  are  calcu¬ 
lated.  These  data  can  be  used  as  inputs  for  a  finite -element  analy¬ 
sis  of  the  fuselage  section. 

14.  Atluri,  S.,  PETROS  3:  A  FINITE-DIFFERENCE  METHOD  AND  PROGRAM  FOR  THE 
CALCULATION  OF  LARGE  ELASTIC-PLASTIC  DYNAMICALLY- INDUCED  DEFORMATIONS 
OF  MULTILAYER  VARIABLE  THICKNESS  SHELLS,  BRL,  U.S.  Army  Aberdeen  Re¬ 
search  &  Development  Center,  Aberdeen  Proving  Ground,  Maryland,  Con¬ 
tract  #DAAD05-68-C-03l4,  Nov.  1971. 

The  governing  equations  for  the  arbitrarily  large-deformation  elastic 
plastic  transient  responses  of  variable-thickness,  hard-bonded,  multi 
layer,  multimaterial,  thin,  Kirchhoff  shells  of  any  initial  shape  are 
formulated  and  solved  by  the  finite-difference  technique.  The 
material  is  assumed  to  be  initially  isotropic  and  to  exhibit  elastic, 
strain-hardening,  strain-rate-sensitive,  and  temperature- dependent 
behavior.  The  structure  may  be  subjected  to  a  variety  of  initial 
velocity  distributions,  transient  mechanical  loads,  and/or  transient 
thermal  loads.  These  capabilities  and  features  are  contained  in  a 
computer  program,  PETROS  3,  which  has  been  applied  to  a  variety  of 
example  problems. 

Included  Is  a  FORTRAN  IV  listing  and  a  description  of  PETROS  3  to¬ 
gether  with  the  data  input  and  solution  output  for  several  example 
problems. 

15.  Haftka,  R.T.,  A  KOITER-TYFE  METHOD  FOR  FINITE  ELEMENT  ANALYSIS  OF 
NONLINEAR  STRUCTURAL  BEHAVIOR,  AFFDL-TR-70-130,  Vol.  I,  Air  Force 
Flight  Dynamics  Laboratory,  Wright-Patterson  Air  Force  Base, 

Ohio,  Nov.  1970. 

Koiter's  method  for  the  asymptotic  analysis  of  post-buckling  behavior 
is  reformulated  in  finite  element  notation  for  application  to  struc¬ 
tures  idealized  by  finite  element  models.  Originally  restricted  to 
the  analysis  either  of  structures  exhibiting  bifurcation  buckling  or 
of  slightly!  imperfect  versions  of  such  structures,  Koiter’s  method 
is  therein  adapted  to  a  more  general  class  of  structures  exhibiting 
the  more  common  snap-through  (limit  point)  type  of  buckling.  This 
adaption  of  Koiter’s  method  is  referred  to  as  the  Modified  Structure 
method.  It  is  accomplished  by  modification  of  the  actual  energy 
functional  to  create  a  hypothetical  modified  structure  having  a 
strictly  linear  pre-buckling  path  along  which  buckling  must  be  of 
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the  bifurcation  type.  Hie  analysis  of  the  actual  structure  is  then 
accomplished  by  application  of  Koiter's  method  through  consideration 
of  the  actual  structure  as  an  imperfect  version  of  the  modified 
structure . 

In  this  way,  the  Modified  Structure  method  operates  within  the 
theoretical  framework  established  by  Koiter,  and  the  effects  of 
prebuckling  nonlinearity  are  approximated  asymptotically.  Various 
levels  of  approximation  are  considered.  Additionally,  the  use  of 
the  Modified  Structure  method,  in  conjunction  with  direct  methods  of 
nonlinear  analysis,  is  examined.  A  highly  accurate  finite  element 
representation  is  employed  in  presenting  a  comprehensive  numerical 
evaluation  of  uhe  Modified  Structure  method  of  analysis  on  the  basis 
of  a  number  of  planar  frame  problems.  Collectively,  these  examples 
exhibit  a  broad  spectrum  of  nonlinear  behavior  characteristics. 
Emphasis  throughout  is  placed  upon  assessing  the  limitations  and  at¬ 
tributes  of  the  Modified  Structure  method  of  analysis.  Conclusions 
regarding  applicability  and  performance  emerge  from  detailed  examina¬ 
tion  of  the  results  obtained. 

1 6.  Stllwell,  W.C.,  and  Ball,  R.E.,  A  DIGITAL  COMPUTER  STUDY  OF  THE 
BUCKLING  OF  SHALLOW  SPHERICAL  CAPS  AND  TRUNCATED  HEMISPHERES,  NASA 
CR  1998,  June  1972. 

Several  user-oriented  digital  computer  programs  for  the  static  analy¬ 
sis  of  shells  of  revolution  exist.  A  detailed  discussion  of  most  of 
these  programs  is  given  in  Reference  1  of  the  report.  Of  particular 
interest  here  is  the  program  developed  by  Ball  (Reference  2)  for  the 
geometrically  nonlinear  analysis  of  arbitrarily  loaded  shells  of  re¬ 
volution.  This  program  is  an  equilibrium  program;  that  is,  it  solves 
for  the  displacement  and  stress  resultant  fields  for  an  arbitrary 
loading  condition.  Since  geometric  nonlinearities  are  included,  the 
magnitude  of  load  that  leads  to  a  condition  of  instaoility  can  be 
determined. 

The  utility  of  the  program  would  be  considerably  enhanced  if  it 
could  be  used  to  determine  bifurcation  buckling  loads  and  the 
behavior  of  the  shell  in  the  vicinity  of  the  bifurcation  load.  This 
latter  feature  is  often  referred  to  as  the  imperfection  sensitivity 
of  the  shell  to  the  load.  As  a  consequence,  the  objective  of  this 
study  was  to  use  the  computer  program  to  examine  the  buckling 
behavior  of  several  shells  subjected  to  axisymmetric  and  nearly 
axisymmetric  loads.  It  was  anticipated  that  an  examination  of  the 
effects  of  the  small  asymmetric  perturbations  upon  the  stability  of 
the  shell  would  disclose  the  bifurcation  buckling  load  and  provide 
a  quantitative  evaluation  of  the  imperfection  sensitivity  of  the 
shell  to  the  load. 

17.  Witmer,  E.A. ,  LARGE  DYNAMIC  DEFORMATIONS  OF  BEAMS,  RINGS,  PLATES  AND 
SHELLS,  AIAA  Journal,  Vol.  I,  No.  2,  August  1963. 
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The  axiaymmetric  responses  of  shells ,  plates,  rings,  and  beams  to 
Impulsive  or  blast  loading  that  produces  large  deformations  involving 
both  the  elastic  and  plastic  regions  of  material  behavior  are  analys¬ 
ed.  A  general  numerical  method  that  includes  (l)  elastic,  (2)  per¬ 
fectly  plastic,  (3)  elastic,  strain -hardening,  and/or  (4)  elastic, 
strain -hardening,  strain-rate  sensitive  material  behavior  and  large 
structural  deflections  is  formulated  and  applied.  In  the  timewise 
sbo-by-step  numerical  analysis,  the  increments  in  stress  resultants 
and  stress  couples  are  determined  by  idealizing  the  shell  thickness 
as  consisting  of  n  concentrated  layers  of  materials  separated  by  a 
material  that  cannot  carry  normal  stresses  but  has  infinite  shear 
rigidity.  The  influences  of  the  number  of  layers  employed  in  the 
idealized  model,  as  well  as  the  aforementioned  various  types  of 
material  behavior,  are  demonstrated.  Theoretical  predictions  of 
time-history  responses  and/or  final  structural  deformations  are  com¬ 
pared  with  experimental  data  for  impact-loaded  spherical  shells,  for 
blast-loaded  circular  plates,  and  for  explosively  loaded  circular 
rings  and  clamped  beams. 

18.  Perrone,  N.,  OiJ  A  SIMPLIFIED  METHOD  FOR  SOLVING  IMPULSIVELY  LOADED 
STRUCTURES  OF  RATE- SENSITIVE  MATERIALS,  Office  of  Naval  Research , 
Washington,  D.C.,  Journal  of  Applied  Mechanics,  A.  S.M.E.,  September  1965 . 

In  an  attempt  to  assess  more  completely  rate-sensitivie  material 
effects,  two  fundamental  structural  elements  are  analyzed:  a  wire 
with  an  impulsively  loaded  end  mass,  and  an  impulsively  loaded  ring. 
The  ring  and  wire  are  made  of  perfectly  plastic,  rate-sensitive 
materials.  In  each  case,  exact  and  approximate  solutions  are  obtain¬ 
ed  for  an  exponential  rate-sensitivity  law.  The  results  suggest 
that  very  good  approximations  to  the  exact  solutions  nay  be  found  by 
utilizing  a  rate-insensitive  material  with  constant  yield  stress 
equal  to  the  initial  dynamic  yield  stress. 

19.  Bodner,  S.R.,  and  Symonds,  P.S.,  EXPERIMENTAL  AND  THEORETICAL  IN¬ 
VESTIGATION  OF  THE  PLASTIC  DEFORMATION  OF  CANTILEVER  BEAMS  SUBJECTED 
TO  IMPULSIVE  LOADING,  Brown  University,  Journal  of  Applied  Mechanics, 
December  1962. 

The  experimental  techniques  and  the  results  obtained  in  a  program  to 
evaluate  the  assumptions  of  dynamic,  rigid-plastic  theory  of  beams 
are  presented.  The  experiments  use  steel  and  aluminum- alloy  canti¬ 
lever  beams  subjected  to  either  a  rapid  velocity  change  at  the  base 
or  to  an  impulsive  load  at  the  tip.  A  rigid-plastic  theory  that  in¬ 
cludes  the  strain-rate  dependence  of  the  yield  stress  and  geometry 
changes  is  outlined  for  the  case  of  tne  tip  impulsive  loading.  Pre¬ 
dictions  of  this  theory  are  in  satisfactory  agreement  with  the  exper¬ 
imental  results. 


10 


tmm mm 


A 


4 


r 


20.  Hibbit ,  H.D.,  et  al,  A  FINITE  ELEMENT  FORUMLATION  FOR  PROBLEMS  OF 
LARGE  STRAIN  AND  LARGE  DISPLACEMENT,  Brown  University,  Int'l  Journal 
of  Solids  and  Structures ,  1970,  Vol  6,  pp.  1069-1086. 

An  incremental  and  piecewise  linear  finite  element  theory  is  developed 
for  the  large-displacement,  large-strain  regime  with  particular  refer¬ 
ence  to  elastic-plastic  behavior  in  metals.  The  resulting  equations, 
though  more  complex,  are  in  a  similar  form  to  those  previously  devel¬ 
oped  for  large-displacement,  small-strain  problems,  the  only  addition¬ 
al  term  being  an  initial  load  stiffness  matrix  which  is  dependent 
on  current  loads.  This  similarity  in  form  means  that  existing  non¬ 
linear  general-purpose  programs  may  easily  be  extended  to  include 
finite  strains.  A  large-displacement,  small-strain  formulation  (as 
applicable  to  problems  of  structural  stability)  is  obtained  from  this 
theory  by  assuming  that  changes  in  length  of  line  elements  and  rela¬ 
tive  rotation  of  orthogonal  line  elements  are  negligible  compared  to 
unity.  The  simplified  equations  are  in  essential  agreement  with  pre¬ 
vious  formulations  in  the  literature.  The  only  difference  which  is 
observed  is  the  persistence  of  the  initial  load  stiffness  matrix, 
which  may  be  significant  in  some  cases. 

21.  Toridio,  T.G.,  and  Khozeimeh,  K.,  INELASTIC  RESPONSE  OF  FRAMES  TO 
DYNAMIC  LOADS,  Journal  of  Engineering  Mechanics,  June  1971* 

A  procedure  of  analysis  is  presented  for  determining  the  elastic- 
inelastic  responsf  of  framed  structures  under  dynamic  loads.  Appli¬ 
cation  of  Hamilton's  principle  in  conjunction  with  the  finite -element 
method  leads  to  the  basic  dynamic  equation  of  the  system  incorporating 
the  plastic  effects  in  the  form  of  equivalent  nodal  forces.  This  ap¬ 
proach  also  allows  the  more  accurate  treatment  of  the  distributed  mass 
of  the  element  than  the  usual  geometrical  method  of  lumping. 

22.  McDaniel,  T.J.,  DYNAMICS  OF  STIFFENED  CYLINDRICAL  SHELLS  WITH  SPATI¬ 
ALLY  VARYING  CURVATURE,  University  of  Dayton  Research  Institute,  Air 
Force  Materials  Laboratory  Report  AFML-TR-72-13^ ,  Air  Force  Systems 
Command,  Wright-Patterson  Air  Force  Base,  Ohio,  July  1972. 

This  report  presents  the  results  of  a  theoretical  study  of  the  effects 
of  spatially  vrrying  curvature  on  the  dynamics  of  cylindrical  panels 
and  stiffened  cylindrical  shell  structures.  The  first  section  of  the 
report  contains  a  brief  discussion  of  the  general  problem  area.  Fol¬ 
lowing  this  discussion,  the  analytical  techniques  for  solving  constant- 
curvature  cylindrical  panels  and  stiffened  cylindrical  shell  problems 
by  a  transfer  matrix  approach  are  reviewed.  These  techniques  are 
found  to  apply  directly  to  the  varying  curvature  shell  analysis,  pro¬ 
vided  the  transfer  matrix  for  this  type  of  shell  can  be  obtained.  An 
analytical  approach  to  obtaining  the  transfer  matrix  for  a  shell  with 
varying  curvature  is  explored.  A  solution  to  the  transfer  matrix  for  a 
cylindrical  shell  with  exponentially  varying  curvature  is  obtained.  In 
a  following  section,  the  preparation  of  this  solution  to  obtain  numeri¬ 
cal  transfer  matrix  is  discussed.  Several  approximate  and  numerical 
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procedures  for  obtaining  a  transfer  matrix  are  explored.  Finally, 
the  dynamic  responses  of  both  single  panels  and  stringer  stiffened 
cylindrical  structures  with  increasing  and  decreasing  curvature  are 
compared  to  similar  structures  with  constant  curvature. 

23.  Bendix  Corp. ,  Final  Engineering  Report,  ENERGY  ABSORBING  CHARACTER¬ 
ISTICS  OF  CRUSHABLE  ALUMINUM  STRUCTURES  IN  A  SPACE  ENVIRONMENT,  NAS  1- 
CR-65096,  July  1965. 

The  research  effort  detailed  in  this  report  involves  the  basic  objec¬ 
tive  of  obtaining  quantitative  design  data  concerning  the  character¬ 
istics  of  aluminum  honeycomb  materials  when  used  in  high  L/D  ratio, 
crushable,  energy-absorbing  capsules.  Nine  configurations  of  honey¬ 
comb  energy-absorbing  capsules  using  alloy  505 6  are  evaluated.  The 
nine  basic  configurations  incorporate  three  cross-sectional  shapes  of 
high,  medium,  and  low  crush  strength,  each  of  which  was  fabricated 
with  cell  axes  oriented  at  angles  of  0/  15  and  30  degrees  to  the 
capsule  longitudinal  axis.  The  characteristics  which  are  studied 
included  specific  energy,  load  onset  rate,  and  rebound.  Variations 
of  these  characteristics  are  investigated  under  controlled  environ¬ 
mental  conditions. 

The  capsules  were  subjected  to  both  static  and  dynamic  loads,  impact 
velocities  from  5  thru  20  feet  per  second,  and  impact  weights  varying 
from  760  thru  3750  pounds.  The  environmental  extremes  under  which 
the  specimens  were  tested  spanned  the  temperature  range  from  -26o°F 
thru  room  temperature  up  to  +300°F,  and  a  vacuum  of  3  x  10“ 1  TORR. 

24 .  Komhauser,  M.,  STRUCTURAL  EFFECTS  OF  IMPACT,  Spartan  Books,  Inc., 
Baltimore,  Md.,  1964. 

Mechanical  impact  has  been  treated  traditionally,  in  the  United 
States,  as  a  subject  which  is  variously  appended  to  college  textbooks 
on  elasticity,  strength  of  materials,  engineering  mechanics,  or 
vibrations.  The  elasticity  and  strength  of  materials  texts  are 
generally  concerned  with  local  surface  effects  or  with  stress  waves; 
the  mechanics  texts  usually  treat  the  whole-body  motions  on  impact; 
and  the  vibrations  texts  are  apt  to  consider  impact  a  special  case 
of  unsteady  vibrations  and  amplification  factors.  In  the  actual  case 
of  impact,  all  of  these  phenomena  and  effects  come  into  play,  and 
the  significance  of  each  effect  must  be  emphasized  relative  to  the 
purpose  of  the  analysis. 

In  this  book  the  emphasis  is  placed  on  go  or  no-go  behavior,  surviv¬ 
ability,  or  failure.  Loading  and  response  must,  or  course,  be  ana¬ 
lyzed.  Nevertheless,  wherever  possible,  the  object  is  to  permit 
estimates  of  failure  directly  in  terms  of  the  loading  conditions. 
Impact  is  a  complex  process.  Given  the  loading  history  of  a  struc¬ 
ture  (and  assuming  no  interaction  with  the  structure’s  response),  it 
is,  in  general,  impossible  to  trace  the  stress  waves  and  their 
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reflections  throughout  the  structure,  resulting  in  vibrations  and 
permanent  set  or  failure.  Many  solutions  of  loading  and  response, 
as  well  as  theories  of  failure  in  terms  of  material  properties,  are 
available  in  the  literature  for  various  idealized  conditions  and  con¬ 
figurations.  What  is  desired  by  the  practicing  engineer  are  some  re¬ 
latively  simple  approaches  to  prediction  of  failure.  This  book  is 
intended  as  a  start  in  this  direction. 

In  attempting  to  provide  engineering  answers  to  some  very  complex 
problems,  liberal  use  has  neen  made  of  "engineering  judgment"  and 
idealizations.  The  soundness  of  each  assumption  is  in  proportion 
to  the  number  of  exact  theoretical  and  experimental  results  available 
for  very  similar  situations.  For  this  reason  it  is  anticipated  that 
modifications  of  the  approaches  presented  therein  will  be  appropriate 
as  the  various  disciplines  produce  more  data. 

The  book  is  organized  in  three  sections:  "Loading  Conditions",  "Re¬ 
sponse  and  Failure  of  Structures",  and  "Effects  of  Impact."  The  sec¬ 
tion  on  response  and  failure  is  expository  in  nature,  introducing 
the  background  theory  and  application  of  the  sensitivity  method  of 
presenting  and  predicting  inertial  failure,  as  well  as  discussing  low- 
speed  and  hypervelocity  impact  effects.  Having  accepted  the  ap¬ 
proaches  recommended  in  the  second  section,  the  practicing  engineer 
may  use  the  curves  and  tables  of  the  first  and  third  sections  for  pre¬ 
diction  of  impact  effects  in  terms  of  the  environmental  input  functions. 

25.  Fisher,  L,  J.,  Jr.,  IANDING- IMPACT-DISSIPATION  SYSTEMS,  NASA,  Techni¬ 
cal  Note  D-975;,  December  1961. 

Analytical  and  experimental  investigations  are  made  to  determine  the 
landing-energy-dissipation  characteristics  for  several  types  of 
earth-landing- impact  systems  having  application  to  reentry  vehicles. 

The  areas  of  study  are  divided  into  three  velocity  regions:  (l)  those 
having  primarily  vertical  velocity,  (2)  those  having  both  moderate 
horizontal  and  moderate  vertical  velocity,  and  (3)  those  having  pri¬ 
marily  horizontal  velocity.  The  impact  systems  discussed  are  braking 
rockets,  gas-filled  bags,  frangible  metal  tubing,  aluminum  hone;,  comb, 
balsa  wood,  strain  straps,  and  both  skid  and  skid-rocker  landings  on 
hard-surface  runways  and  on  water. 

The  report  states  that  it  appears  feasible  to  evaluate  landing-gear 
systems  for  reentry  vehicles  by  computational  methods  and  free-body 
landing  techniques  with  energy  dissipation  for  an  earth  landing  of 
such  a  vehicle.  Some  systems  are  more  efficient  thm  others,  some 
package  better  than  others,  and  a  variety  of  promising  systems  are 
under  study.  Horizontal  energy  dissipation  is  simpler  to  deal  with 
than  vertical  energy  dissipation  since  translational  friction  is  all 
that  is  involved;  however,  runout  behavior  becomes  a  factor.  Verti¬ 
cal  velocity  can  also  be  a  big  factor  when  high  flight-path  angles  are 
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associated  with  even  moderate  horizontal  velocities.  High-speed  land¬ 
ings  are  particularly  a  problem,  especially  high-speed  water  landings, 
and  indications  are  that  if  large  horizontal  velocities  are  involved  in 
hard-surface  landings,  a  selected  site  will  be  required. 

26.  McGehee,  J.  R.,  A  PRELIMINARY  EXPERIMENTAL  INVESTIGATION  OF  AN  ENERGY- 
ABSORPTION  PROCESS  EMPLOYING  FRANGIBLE  METAL  TUBING,  NASA  Technical 
Note  D-1477,  1962. 

A  highly  efficient  energy-absorption  process,  employing  frangible 
meted  tubing  as  the  working  element,  is  investigated.  A  preliminary 
experimented  investigation  is  conducted  to  determine  the  variation  of 
the  average  fragmenting  stress  of  2024-T3  aluminum- alloy  tubing  with 
the  pertinent  parameters  of  this  process.  Limited  tests  were  made  to 
determine  the  feasibility  of  employing  this  process  in  a  landing-gear 
system.  A  l/ 5-scale  model  of  a  proposed  manned  spacecraft  with  a 
landing  gear  incorporating  this  process  is  employed  in  these  tests. 

The  results  of  this  investigation  show  that  the  fragmenting  process 
produces  a  fluctuating  force  with  displacement,  but  for  a  fixed  set 
of  parameters,  the  force  about  which  the  fluctuation  occurs  is  ap¬ 
proximately  constant.  A  large  force  which  occurs  when  the  process  is 
started  with  the  unaltered  tube  seated  symmetrically  in  the  die  can 
be  reduced  most  effectively  by  tapering  the  wall  thickness  over  a 
short  length  at  the  die  end  of  the  tube.  The  average  fragmenting 
stress,  for  2024-T3  aluminum- alloy  tubing  and  the  range  of  parameters 
investigated,  appears  to  be  independent  of  the  ratio  of  wall  thickness 
to  tube  diameter,  but  varies  as  the  cube  of  the  ratio  of  the  wall 
thickness  to  the  radius  ol  .he  forming  die.  The  fragmenting  stress 
obtained  at  12,000  inches  per  minute  was  about  60  percent  higher  than 
those  obtained  at  1  inch  per  minute.  The  2021* -T3  aluminum- alloy  tub¬ 
ing,  when  fragmenting  on  a  die  at  90  percent  of  the  yield  stress,  is 
capable  of  absorbing  31,000  foot-pounds  of  energy  per  pound  of 
material.  This  energy- absorption  capability  is  greater  than  that  of 
the  most  frequently  considered  processes;  for  example,  the  crushing 
of  balsa  wood,  aluminum  honeycomb,  or  pressurized  thin-walled  metallic 
cylinders.  Model  tests,  employing  frangible  tubing  as  the  working 
element  in  the  landing  gear,  indicate  that  this  process  is  suitable 
for  use  in  a  lead-alleviation  application. 

27.  Kroell,  C.  K.,  A  SIMPLE,  EFFICIENT,  ONE  SHOT  ENERGY  ABSORBER,  General 
Motors  Research  Laboratory,  Warren,  Michigan,  Shock,  Vibration  and 
Associated  Environments,  Part  III,  Bulletin  No.  30,  19&2. 

This  paper  describes  a  single-shot  expendable  energy  absorber  which 
has  been  developed  at  the  General  Motors  Research  Laboratory.  The 
device  is  inherently  simple  and  is  characterized  by  a  rectangular 
force-displacement  relationship  and  high  specific  energy  absorption 
capacity.  Both  a  qualitative  discussion  of  the  mechanics  of  the 


plastic  deformation  process  involved  and  a  graphical  summary  of  the 
experimental  performance  data  which  have  been  collected  to  date  are 
presented. 


29. 


Weinberg,  L.W.T.,  and  Tumbov,  J.W.,  Fh.D.,  SURVIVABILITY  SEAT  DESIGN 
DYNAMIC  TEST  PROGRAM,  Aviation  Safety  Engineering  and  Research, 
USAAVLABS  Technical  Report  65-^3 t  U.S.  Army  Aviation  Materiel  Labora¬ 
tories,  Fort  Eustis,  Virginia,  July  19&5»  AD  621718. 


This  report  presents  the  results  of  a  series  of  dynamic  tests  conduct¬ 
ed  with  four  different  concepts  of  experimental  ere-  seats. 


The  experimental  seats  were  designed  and  constructed  by  four  heli¬ 
copter  manufacturers.  The  seats  were  designed  to  withstand  static 
load  factors  equivalent  to  those  recommended  in  TRECOM  Technical  Re¬ 
port  63-^,  "Crew  Seat  Design  Criteria  for  Army  Aircraft",  dated 
February  1963. 


The  design  load  factors  recommended  in  the  above-referenced  report 
are  as  follows:  longitudinal  -  45G  for  0.10  second;  lateral  -  L5G 
for  0.10  second;  and  vertical  -  25G  for  0.10  second. 


Special  kits  for  small-arms  ballistic  protection  were  also  designed 
and  installed  in  the  seats  tested. 


These  seats  were  designed  exclusively  using  static  load  factors.  No 
previous  testing  was  conducted  by  any  seat  manufacturer  prior  to  the 
conduct  of  these  tests. 


The  four  seats  were  tested  under  four  load  conditions.  Two  of  the 
conditions  involved  simultaneous  half  loads  on  the  seats  in  two  dif¬ 
ferent  seat  positions,  and  two  of  the  conditions  involved  full  loads 
in  two  different  seat  positions. 


Only  one  of  the  four  seats  tested  withstood  the  loads  imposed  for 
all  four  conditions.  Three  of  the  seats  failed  and  were  damaged 
beyond  economical  repair  when  each  was  subjected  to  the  first  full¬ 
load  test  condition. 


This  report  also  includes  a  detailed  description  of  an  acceleration 
device  which  was  specifically  designed  and  fabricated  for  this  series 
of  tests. 


Langhaar,  H.  L.,  THEORETICAL  AND  EXPERIMENTAL  INVESTIGATIONS  OF  THIN- 
WEBBED  PLATE-GIRDER  BEAMS,  Transactions  of  the  ASME .  October  I9U3. 


A  simple,  semirational  theory  for  the  design  of  webs  and  flange 
rivets  of  thin-webbed  rectangular  plate-girder  shear  beams  is  pie- 
sen  ted.  Calculations  of  shear  loads  to  cause  web  rupture  and  flange 
rivet  failure  are  compared  with  test  data  from  ,27  beams. 
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30.  Perry,  D.  J.,  AIRCRAFT  STRUCTURES,  McGraw  Hill  Book  Co.,  New  York, 

1950. 

In  this  hook  an  attempt  is  made  to  emphasize  basic  structural  theory 
which  will  not  change  as  new  materials  and  new  construction  methods 
are  developed.  Most  of  the  theory  is  applicable  for  any  design  re¬ 
quirements  and  for  any  materials.  The  design  engineer  may  then  sup¬ 
plement  this  theory  with  the  detail  design  specifications  and  the 
material  properties  which  are  applicable  to  his  particular  airplane. 

Heavy  emphasis  is  placed  on  the  application  of  the  elementary  prin¬ 
ciples  of  mechanics  to  the  analysis  of  aircraft  structures. 

31.  Jones,  N.,  et  al,  THE  DYNAMIC  PLASTIC  BEHAVIOR  OF  SHELLS,  MIT  Report 
71-6,  1971. 

A  survey  is  made  of  the  literature  published  previously  on  the  in¬ 
elastic  behavior  of  shells  subjected  to  dynamic  loads.  An  experi¬ 
mental  investigation  is  also  being  undertaken  to  examine  the  behavior 
of  various  cylindrical  shell  panels  which  are  loaded  with  an  impulse 
on  the  inner  surface.  The  panels  are  fully  clamped  along  the  two 
longitudinal  edges  and  free  on  the  other  two.  The  initial  kinetic 
energy  of  the  dynamic  loads  is  sufficiently  large  to  cause  inelastic 
behavior  and  to  produce  maximum  permanent  transverse  deflections  of 
up  to  nearly  twice  the  corresponding  panel  thickness.  Tests  are  con¬ 
ducted  cn  mild  steel  and  aluminum  606I-T6  panels  which  have  various 
thicknesses  and  included  angle  of  90°  approximately. 

32.  Mitchell,  B.,  THE  DYNASORB  ENERGY  ABSORBER,  Lockheed  Report  LR  16735, 
March  19^3. 

This  report  introduces  a  method  of  absorbing  energy  while  maintaining 
a  constant  load  level.  The  tubular  load-carrying  member  can  support 
its  full  load  capabilit1,  and  still  be  almost  totally  consumed  by  roll¬ 
ing  up  one  end.  The  energy  curves  presented  are  typical  and  show  ex¬ 
cellent  load  consistency  with  efficiencies  as  high  as  450,000  in.- 
lb/lb. 

33.  Mitchell,  B.,  SHOCK  ABSORPTION  WITH  ONE  SLOT  TUBES,  Lockheed  Report 
LR  16369,  June  1963. 

The  work  covered  in  this  report  is  a  follow-on  of  an  independent  re¬ 
search  program  that  developed  a  very  efficient  method  of  absorbing 
energy  by  load  control  and  end  roll-up  of  tubes.  This  study  investi¬ 
gates  12  different  materials  and  improves  the  control  rings  to  give  a 
smooth,  nearly  rectangular,  energy  curve.  The  load  level  can  be  set, 
by  suitable  design,  at  any  desired  point  within  the  column  or  local 


strength  limits  of  the  tube.  The  maximum  efficiency  attained  here 
with  a  2-inch-diameter  by  .049  wall  4130  steel  tube  heat  treated  to 
200,000  psi  is  600,000  in. -lb  per  lb  of  tube  weight.  This  represents 
a  mean  compressive  stress  of  189,000  psi  over  the  total  6  inches  of 
travel, 

3‘i.  Mitchell,  B.,  DESIGN  NOTES  FOR  THE  DYNASORB  ENERGY  ABSORBER,  Lockheed 
Report'  LR  17201,  December  1963* 

The  Dynasorb  Energy  Uhit,  developed  in  the  Lockheed  Engineering 
Laboratory,  is  adaptable  to  many  applications  as  a  shock  absorber  or 
load  limiter.  Desigi  procedures  based  on  previously  reported  data 
are  described,  and  several  typical  installations  are  illustrated. 

This  report  is  submitted  in  fulfillment  of  the  reporting  requirements 
of  a  1963  Independent  Development  Project,  "Energy  Absorption  Pro¬ 
ducts  . " 

35.  Mitchell,  B.,  ENERGY  ABSORPTION  AT  HIGH  SPEED  VERTICAL  LANDING, 
Lockheed  Report  LR  21023,  November  1967* 

This  report  describes  the  use  of  the  Lockheed  Dynasorb  energy  absorb¬ 
er  in  landing  impact  tests  of  a  3400-lb  simulated  vehicle  structure. 
The  structure,  released  from  a  helicopter  and  tethered  by  parachute, 
impacted  at  a  velocitv  of  112  feet  per  second. 

36.  Per rone,  N.,  RESPONSE  OF  RATE  SENSITIVE  FRAMES  TO  IMPULSIVE  LOAD, 

ASME  Journal  of  Applied  Mechanics,  February  1971 ,  pp.  49-62. 

A  relatively  convenient  method  is  presented  herein  to  determine  the 
plastic  response  of  rate-sensitive  frameworks.  The  technique  repre¬ 
sents  a  significant  extension  of  usual  limit  analysis  type  approaches 
as  well  as  similar  efforts  applied  to  rate  insensitive  structures. 


The  method  consists  of  assuming  a  modal  deformation  pattern  for  the 
structure  in  question,  determination  of  the  initial  modal  velocity 
utilizing  the  appropriate  criteria,  and  integration  in  some  form  of 
the  equations  of  motion  after  estimating  the  magnitude  of  the  dynamic 
yield  moments  at  the  rate-sensitive  plastic  hinge  locations.  It  is 
recommended  that  an  effective  hinge  length  of  the  order  of  three  beam 
depths  should  suffice  in  estimating  the  strain  rate  magnitude.  For 
the  perfectly  plastic  case,  the  dynamic  yield  moment  is  assumed  to  be 
a  constant  with  time. 

A  square  portal  frame  made  of  a  solid  rectangular  cross-section  and 
loaded  under  a  horizontal  impulsive  load  is  considered  in  some  detail. 
It  is  shown  that  a  critical  initial  velocity  exists. 

The  effects  of  strain  hardening  and  pulse  load  application  are  con¬ 
sidered.  Minor  modifications  of  the  impulsively  loaded,  perfectly 
plastic  situation  are  necessary  to  accommodate  these  extensions. 
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The  usual  limitations  should  be  noted.  It  is  assumed  that  buckling 
will  not  occur,  that  the  median  surface  or  membrane  effects  are 
negligible,  and  that  elastic  effects  are  negligibly  small  relative  to 
plastic  flow  because  the  former  are  omitted. 

Portal  frame  experiments  under  high  intensity  loading  similar  to 
those  conducted  at  Brown  University  on  cantilever  beams  would  be 
most  welcome.  Most  of  the  experiments  performed  to  date  have  been 
in  a  much  lower  load  range  where  elastic  and  plastic  effects  are  of 
comparable  order.  The  true  limiting  strength  of  the  framework  could 
be  more  completely  assessed  only  if  larger  loads  are  applied.  Other 
tests  have  been  reported,  and  hopefully  these  results  will  soon  be 
available . 

37.  Jones,  N.,  INFLUENCE  OF  STRAIN-HARDENING  AND  STRAIN-RATE  SENSITIVITY 
ON  THE  PERMANENT  DEFORMATION  OF  IMPULSIVELY  LOADED  RIGID  PLASTIC 
BEAMS,  International  Journal  of  Mechanical  Science,  1967,  Vol.  9, 
pp.  777-79^ 

A  simple  method  is  presented  for  estimating  the  combined  influence  of 
strain-hardening  and  strain-rate  sensitivity  on  the  permanent  defor¬ 
mation  of  rigid-plastic  structures  loaded  dynamically.  A  study  is 
made  of  the  particular  case  of  a  beam  supported  at  the  ends  by  im¬ 
movable  frictionless  pins  and  loaded  with  a  uniform  impulse.  The 
results  of  this  work  indicate  that,  when  stress -hardening  or  strain- 
rate  sensitivity  are  considered,  permanent  deformations  are  experi¬ 
enced  which  are  similar  to  those  predicted  by  an  analysis  retaining 
both  effects  simultaneously. 

38.  Ni,  C.  M. ,  IMPACT  RESPONSE  OF  CURVED  BOX  BEAM-COLUMNS  WITH  LARGE 
GLOBAL  AND  LOCAL  DEFORMATIONS,  General  Motors  Research  Laboratory, 
Warren,  Michigan,  AIAA  Paper  73-^01. 

A  numerical  approach  based  on  a  lumped-mass  model  is  developed  for 
investigating  the  impact  response  of  curved  box  beam-columns  with 
large  global  and  local  deformations.  An  empirical  formula  which  re¬ 
lates  the  changes  of  depth  and  bending  angle  of  a  beam  cross-section 
is  obtained  to  take  into  account  the  local  deformations  of  the  cross- 
sections.  The  strain-hardening  and  the  strain-rate  properties  of  the 
material  are  considered  in  this  analysis.  The  correlation  between 
the  present  analysis  and  test  results  is  very  good.  The  results 
obtained  indicate  that  the  strain -wave  propagations  due  to  the  impact 
and  the  strain-rate  sensitivity  of  the  material  play  the  key  roles 
In  increasing  the  energy-absorbing  capacity  of  the  structure  when 
subjected  to  high-speed  impact. 


18 


'  1  ■v" '»»'»■  *»"«  .< 


39. 


4o. 


■bmwwbwww  v  i  n  111— .  wmrnmmimmmmmssJmMmMMmammmmaaxr — JR9HT 

O’Bryan,.  T.C.,  and  Hatch,  H.G. ,  Jr.,  LIMITED  INVESTIGATION  OF  CRUSH- 
ABLE  STRUCTURES  FOR  ACCELERATION  PROTECTION  OF  OCCUPANTS  OF  VEHICLES 
AT  LOW  IMPACT  SPEEDS,  NASA  Technical  Note  D-I58,  1958. 


A  limited  investigation  is  made  to  determine  the  characteristics  of 
three  materials  to  see  hov  they  can  he  applied  for  human  protection 
against  accelerations  encountered  at  low  impact  speeds.  As  a  result, 
if  given  man’s  physiological  tolerance  to  abrupt  acceleration,  which 
has  not  yet  been  well-defined,  an  alleviation  system  can  be  designed. 

Foamed  plastics  require  considerable  depth  to  provide  a  given  stop¬ 
ping  distance  for  impact  alleviation,  and.  their  use  will  require 
some  control  of  rebound.  They  can  be  made  soft  enough  to  obtain  the 
low  onset  of  acceleration  that  may  be  necessary  for  man  where  depth 
is  not  limited. 

Aluminum  honeycomb  is  an  efficient  material  for  impact  load  allevia¬ 
tion  from  the  standpoint  of  usable  material  depth,  and  it  exhibits 
very  little  rebound.  The  stiffness  of  the  material  results  in  a  very 
high  initial  onset  rate  of  acceleration.  For  many  installations  this 
may  be  controlled  by  reducing  the  initial  loading  area  of  contact  to 
get  the  material  to  start  failing. 

Jones,  Norman,  THE  INFLUENCE  OF  LARGE  DEFLECTION  ON  THE  BEHAVIOR  OF 
RIGID-PLASTIC  CYLINDRICAL  SHELLS  LOADED  IMPULSIVELY,  Journal  of 
Applied  Mechanics,  ASME,  June  1970,  pp.  417-425. 


In  order  to  gain  some  insight  into  the  importance  and  influence  of 
finite  deflections  on  the  response  of  shells  loaded  dynamically,  this 
article  studies  theoretically  the  behavior  of  a  cylindrical  shell  sub¬ 
jected  to  a  uniform  axisymmetric  impulsive  pressure.  The  cylindrical 
shell  is  assumed  to  be  made  from  a  rigid,  perfectly  plastic  material, 
and  the  external  energy  imparted  to  the  shell  is  much  greater  than  the 
total  strain-energy  which  can  be  absorbed  elastically.  The  results  of 
the  investigation  indicate  that  geometry  changes  are  important  and 
should  be  retained  when  studying  the  behavior  of  cylindrical  shells 
loaded  dynamically. 

4l.  D’Amato,  R.,  STATIC  POST-FAILURE  STRUCTURAL  CHARACTERISTICS  OF  MULTI¬ 
WEB  BEAMS,  WADC  TR  59-112,  February  1959. 

The  structural  behavior  of  aircraft  components  is  discussed  in  rela¬ 
tionship  to  the  general  lethality  problem,  the  concept  of  post¬ 
failure  restoring  force  for  built-up  structures  of  multiweb  beams  is 
considered,  and  a  relatively  simple  analysis  is  developed  to  compute 
the  static  post-failure  behavior  of  multiweb  beams  in  pure  bending. 

An  experimental  investigation  is  conducted  to  examine  the  validity  of 
the  analysis  and,  within  the  range  of  parameters  considered,  agreement 
between  experiment  and  theory  is  satisfactory.  Both  the  theoretical 
analysis  and  the  experimental  results  indicated  that  an  arbitrary 
extrapolation  of  post-failure  data  on  the  basis  of  ultimate  strength 
can  be  very  misleading. 
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42.  Rawlings,  B.  ENERGY  ABSORPTION  OF  DYNAMICALLY  AND  STATICALLY  TESTED 
MILD  STEEL  BEAMS  UNDER  CONDITIONS  OF  GROSS  DEFORMATION .  International 
Journal  of  Mechanical  Science.  Pergamon  Press. ,  Ltd.,  1967,  Vol.  9, 
pp.  633-649.  Printed  in  Great  Britain. 


An  account  is  given  of  dynamic  and  static  tests  on  mild  steel  members 

deformed  unde:  conditions  of  gross  geometry  change.  An  evaluation  is 

made  of  the  rigid-plastic  theory,  taking  account  of  change  of  geometry,  1 

and  also  the  elastic-plastic  theory,  assuming  deformation  to  occur  in  ^  ; 

the  static  mode,  in  predicting  the  behavior  of  the  members.  *  j 

< 

t 

43.  Thompson,  J.E.,  VEHICLE  CRUSH  PREDICTION  USING  FINITE-ELEMENT  TECH-  j 

NIQUES,  Chrysler  Corp.,  SAE  Paper  73-157.  January  1973* 

The  principal  objective  of  this  investigation  is  to  develop  analyti-  \ 

cal  tools  in  the  form  of  computer  programs  which  will  permit  the  auto-  j 

mobile  designer  to  predict  the  crush  characteristics  of  a  given  car 
structure  due  to  forces  generated  in  a  variety  of  impact  modes.  j 

i 

The  predictive  or  control  capability  is  embodied  in  two  large  computer  j 

programs.  "TELSAP"  forms,  reduces,  and  inverts  the  vehicle  structure  ! 

mass  matrix  expressed  relative  to  a  datum  coordinate  system  and  writes  j 

the  mass  matrix  and  its  inversion  onto  a  file  for  reading  by  the  "CRUSH"  j 

program.  "CRUSH"  is  a  general  matrix  structural  analysis  program  which  1 

calculates  the  large,  plastic,  rate-sensitive  response  of  an  intercon¬ 
nected  beam  structure  due  to  known  dynamic  boundary  displacement  inputs. 

The  theoretical  bases  and  assumptions  employed  in  developing  these  pro¬ 
grams  are  described  along  with  a  detailed  discussion  of  how  the  auto-  { 

mobile  design  engineer  might  use  them  to  develop  a  new  vehicle  j 

structure. 

Experimental  correlation  with  the  computer  models  is  given  for  a 
vehicle-to-vehicle  90-degree  intersection  collision  between  a  special 
rigid  moving  barrier  and  an  intermediate-size  four-door  sedan.  The 
structural  model  is  further  correlated  with  a  laboratory  test  of  a 
clamped-clamped  beam  struck  at  its  center  by  an  impact  pendulum.  The 
correlation  indicates  general  agreement  between  experimental  and  an¬ 
alytical  results. 


By  using  the  computer  programs  developed  in  this  investigation,  the 
automobile  designer  is  able  to  reduce  the  amount  of  testing  required 
to  prove  his  design,  and  is  able  to  identify  the  benefits  of  a  parti¬ 
cular  structural  reinforcement  with  a  minimum  of  development  time  and 
expense . 

44.  Burgmann,  J.  B.,  and  Rawlings,  B.,  DYNAMIC  PLASTIC  ANALYSIS  OF  PIN- 
JOINTED  FRAMES,  International  Journal  of  Mechanical  Science.  Pergamon 
Press,  1968,  Vol.  10,  pp.  967-980,  Printed  in  Great  Britain.  (Revised 
30  July  1968). 
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The  paper  presents  an  analysis  of  a  pin- jointed  frame  subjected  to 
dynamic  or  impulsive  overload,  of  sufficient  magnitude  to  cause  per¬ 
manent  deformation.  Rigid-plastic  behavior  of  tensile  members  is 
assumed  and,  as  a  first  approximation,  a  similar  behavior  of  compres¬ 
sive  members  is  assumed,  although  modifications  to  account  for  other 
characteri sties  are  also  discussed. 

The  behavior  is  considered  in  terms  of  the  kinematic  conditions,  dy¬ 
namic  conditions,  and  the  load- deformation  characteristics  assumed 
for  the  members. 

45.  Martin,  J.B.,  MODE  APPROXIMATION  FOR  IMPULSIVELY  LOADED  STRUCTURES  IN 
THE  INELASTIC  RANGE,  Proceedings  of  the  Southampton  1969  Civil  En¬ 
gineering  Material  Conference. 

A  convergence  approximation  technique,  based  on  the  uniqueness  proof, 
is  reviewed  for  impulsively  loaded  rigid  plastic  and  rigid  visco¬ 
plastic  structures.  Emphasis  is  given  to  the  use  of  mode  or  quasi¬ 
mode  solutions  and  their  usefulness  in  establishing  a  general  approxi¬ 
mating  procedure  and  in  providing  insight  into  the  important  aspects 
of  the  gross  structural  response. 

46.  Jensen,  W.R.,  Plaby,  W.E.,  and  Prince,  N.,  MATRIX  ANALYSIS  METHODS  FOR 
ANISOTROPIC  INELASTIC  STRUCTURES,  AFFDL-TR-65-220,  April  19 66. 

Most  aerospace  structural  materials  exhibit  some  degree  of  aniso¬ 
tropic  strain -hardening.  During  the  past  few  years,  several  methods 
have  appeared  in  the  literature  for  introducing  inelastic  isotropic 
material  behavior  effects  into  existing  matrix  analysis  procedures 
using  the  incremental  theory  of  plasticity.  A  review  is  presented  of 
these  methods  and  a  step-by-step  routine  known  as  the  "Constant 
Strain"  method  is  selected  for  the  development  of  an  anisostropic  in¬ 
elastic  procedure. 

47.  Isaakson,  G. ,  Armen,  H.,  Jr.,  and  Pipko,  A.,  DISCRETE  ELEMENT  METHODS 
FOR  THE  PLASTIC  ANALYSIS  OF  STRUCTURES,  NASA  CR  803,  October  1967. 

This  study  deals  with  the  extension  of  finite-element  methods  to  pro¬ 
vide  analytical  means  for  determining  the  failure  loads  of  aeronauti¬ 
cal  structures.  Two  areas  are  considered  as  related  to  predicting 
failure  loads:  inelastic  stress  analysis  in  the  presence  of  load 
cycling,  and  plastic  buckling  of  the  bifurcation  type. 

Finite-element  inelastic  stress  analysis  methods  are  extended  to 
take  into  account  the  Bauschinger  effect  for  biaxial  stress  states 
using  a  plasticity  theory  based  on  Ziegler’s  modification  to  Prager's 
kinematic  hardening  theory.  This  methodology  is  applied  to  several 
structures  representative  of  aeronautical  construction,  including  a 
notched  plate,  a  shear  lag  specimen,  and  a  swept  wing.  Good  corre¬ 
lation  is  obtained  between  analytical  and  experimental  results  for 


the  strains  at  the  root  of  the  notched  plate  subjected  to  load 
cycling  in  the  plastic  range. 

Finite-element  buckling  methods  are  also  extended  to  consider  plastic 
buckling  using  Stowell ' s  formulation  for  implementing  a  deformation 
plasticity  theory  into  the  buckling  theory.  Sample  calculations  are 
carried  out  for  the  plastic  buckling  of  a  flat  plate  with  various 
geometries  and  edge  conditions. 

48.  Stricklin,  J.A.,  et  al,  NONLINEAR  DYNAMIC  ANALYSIS  OF  SHELLS  OF  REVO¬ 
LUTION  BY  MATRIX  DISPLACEMENT  METHOD,  AIAA  Journal.  Vol.  9,  No.  4, 
April  1971 ,  p.  629. 

A  formulation  and  computer  program  is , developed  for  the  geometrically 
nonlinear  dynamic  analysis  of  shells  of  revolution  under  symmetric 
loads.  The  nonlinear  strain  energy  expression  is  evaluated  using 
linear  functions  for  all  displacements.  Five  different  procedures 
are  examined  for  solving  the  equations  of  equilibrium,  with  Houbolt's 
method  proving  to  be  the  most  suitable.  Solutions  are  presented  for 
the  symmetrical  and  asymmetrical  buckling  of  shallow  caps  under  step 
pressure  loadings  and  a  wide  variety  of  other  problems,  including 
some  highly  nonlinear  ones. 

49.  Owens,  R.  H.,  and  Symonds,  P.S.,  PLASTIC  DEFORMATIONS  OF  A  FREE  RING 
UNDER  CONCENTRATED  DYNAMIC  LOADING,  ASME  Journal  of  Applied  Mechanics  , 
December  1955#  P*  524. 

A  concentrated  time-dependent  force  acts  on  an  unsupported  thin  ring 
along  a  diameter.  The  problem  considered  in  this  paper  is  to  deter¬ 
mine  the  deformations  of  the  ring  when  the  force  magnitudes  are  such 
that  plastic  strains  occur  which  are  large  compared  to  the  elastic 
strains.  By  neglecting  elastic  strains  and  assuming  Ideally  plastic 
behavior,  approximations  to  the  final  deformations  of  the  ring  are 
obtained.  The  analysis  is  developed  for  force  pulses  of  arbitrary 
shape,  but  numerical  results  are  obtained  only  in  the  special  case 
of  a  rectangular  force  pulse.  A  criterion  is  stated  for  conditions 
when  this  type  of  analysis  can  be  expected  to  provide  satisfactory 
results. 


50.  Lee,  H.,  and  Symonds,  P.S.,  LARGE  PLASTIC  DEFORMATIONS  OF  BEAM  UNDER 
TRANSVERSE  IMPACT,  ASME  Journal  of  Applied  Mechanics.  September 
1952,  p.  308. 

A  comparatively  simple  method  of  analysis  is  developed  to  determine 
the  deformations  in  a  beam  subjected  to  lateral  impact  of  such  a 
magnitude  that  plastic  strains  occur  which  are  large  compared  with 
elastic  strains.  A  useful  approximation  to  the  motion  then  can  be 
obtained  by  neglecting  elastic  strains  and  considering  rigid-body 
motion  of  segments  of  the  beam  Joined  at  plastic  hinges  where  the 
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entire  deformation  takes  place.  A  method  of  analyzing  such  a  situa¬ 
tion  is  described  and  applied  to  a  beam  subjected  to  central  impact. 
The  approximate  final  permanent  deformation  is  obtained;  this  includes 
deformation  during  application  of  the  load,  and  plastic  flow  which 
continues  afterward  when  the  kinetic  energy  of  the  motion  generated 
by  the  impact  is  transformed  into  additional  plastic  deformation.  A 
criterion  is  given  for  conditions  when  this  type  of  theory  can  be 
expected  to  provide  a  satisfactory  analysis.  The  method  of  solution 
provides  an  interesting  analogy  to  .the  concept  of  static  determinacy 
which  has  been  used  in  tne  analysis  of  quasi-static  plastic-flow 
problems. 

51.  Prager,  W.,  A  NEW  METHOD  OF  ANALYZING  STRESSES  AND  STRAINS  IN  WORK¬ 
HARDENING  PLASTIC  SOLIDS,  ASME  Journal  of  Applied  Mechanics .  December 
1956,  p.  493. 

For  work-hardening  plastic  solids,  segment-wise  linear  yield  condi¬ 
tions  and  the  associated  flow  rules  constitute  a  reasonable  compro¬ 
mise  between  the  mathematically  convenient  but  physically  unsound 
total  stress-strain  laws  and  the  physically  sound  but  mathematically 
inconvenient  incremental  laws.  They  allow  total  stress-strain  laws  to 
be  used  in  the  small,  but  retain  the  characteristic  features  of  the 
incremental  laws  in  the  large.  The  use  of  a  segmentwise  linear  yield 
condition  and  the  associated  flow  rule  is  illustrated  by  the  analysis 
of  the  bending  moments  and  deflections  of  a  simply  supported  circular 
plate  that  is  made  of  a  work-hardening  material  and  subjected  to  a 
uniformly  distributed  transverse  load. 

52.  Morino,  L.,  Leech,  J.W.,  and  Witmer,  E.A.,  AN  IMPROVED  NUMERICAL  CAL¬ 
CULATION  TECHNIQUE  FOR  LARGE  ELASTIC -PLASTIC  TRANSIENT  DEFORMATIONS 
OF  THIN  SHELLS,  Part  1,  ASME  Journal  of  Applied  Mechanics,  June  1971, 
P.  423. 

In  this  paper,  the  governing  differential  relations  which  describe 
the  large-deflection  elastic-plastic  dynamic  responses  of  arbitrarily 
shaped  thin  Kirchhoff  shells  are  given, including  recent  improvements. 
These  relations  are  then  cast  into  finite-difference  form  for  numeri¬ 
cal  solutions.  These  finite-difference  relations  are  employed  in  a 
computer  program,  PETROS  2,  which  has  been  applied  to  evaluate  the 
"analysis  improvements"  by  comparing  PETROS  2  predictions  with  those 
of  the  earlier  analysis  and  with  experimental  results. 

53.  Pifko,  A.  Issakson,  A  FINITE-ELEMENT  METHOD  FOR  THE  PLASTIC  BUCKLING 
ANALYSIS  OF  PLATES,  Grumman  Aerospace  Corporation,  Bethpage,  N.Y., 

AIAA  Journal  of  Applied  Mechanics.  Vol.  7>  No.  10,  October  19 69. 

An  existing  finite-element  technique  for  elastic  buckling  of  plates 
is  extended  to  include  the  case  of  plastic  buckling.  The  Stowell 
theory  for  the  plastic  buckling  of  flat  plates  is  used  in  conjunction 
with  the  finite-element  technique.  Application  is  made  to  rectangular 
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places,  and  results  are  presented  for  a  variety  of  boundary  support 
conditions  and  several  different  edge  loading  conditions. 

54.  Armen, H. ,  Jr.,  Pifko,  A.,  and  Levine,  H.S.,  FINITE  ELEMENT  ANALYSIS 
OF  STRUCTURES  IN  THE  PLASTIC  RANGE,  NASA  CR  1649,  February  1971* 

The  present  report  is  concerned  with  the  development  of  finite- 
element  methods  for  the  treatment  of  the  nonlinear  behavior  of 
complex  structures.  It  represents  an  extension  of  a  previous  study 
reported  in  NASA  Contractor's  Report  CR-803.  The  nonlinearity  may 
be  of  two  types  -  material  nonlinearity  associated  with  plastic  de¬ 
formation,  and  geometric  nonlinearity  associated  with  the  changing 
geometry  of  the  structure  as  it  deforms  -  or  it  may  involve  a  combin¬ 
ation  of  the  two.  Effects  due  to  creep  and  other  time-dependent 
material  properties  sure  neglected. 

The  methods  developed  are  applicable  to  loading  conditions  that  cause 
membrane  stress  states  or  pure  bending,  or  both  in  combination.  The 
Prsiger-Ziegler  kinematic  hardening  theory  of  plasticity  is  incorporat¬ 
ed  in  the  finite-element  methods  to  allow  for  consideration  of  the 
plastic  response  of  structures  subjected  to  realistic  loading  condi¬ 
tions,  including  cyclic  loadings  that  cause  stress  reversals  into  the 
plastic  range.  Ideally,  plastic  behavior  is  also  included  to  provide 
capability  for  predicting  the  collapse  load  of  structures.  The  plas¬ 
ticity  theory  is  implemented  in  the  finite-element  analysis  by  using 
an  incremental  approach  in  conjunction  with  the  initial  strain  con¬ 
cept,  with  plastic  strains  interpreted  as  initial  strains. 

The  treatment  of  geometric  nonlinearity  requires  use  of  an  increment¬ 
al  technique  in  which  the  internal  forces  and  configuration  of  the 
structure  are  continuously  updated  to  account  for  its  changing  geo¬ 
metry. 

The  methods  developed  are  applied  to  a  number  of  sample  structures. 
For  membrane  stress  states  alone,  the  analysis  employs  a  triangular 
finite  element  in  which  stress  and  strain  vary  linearly.  This  ele¬ 
ment  is  used  for  the  plastic  analysis  of  a  variety  of  structures 
characterized  by  regions  of  rapid  stress  variation  and  subjected  to 
cyclic  loading  resulting  in  reversed  plasticity.  Comparisons  of  the 
results  of  the  analysis  and  experimental  data  indicate  good  correla¬ 
tion. 

Plastic  analyses  are  also  performed  for  a  Variety  of  beam  and  plate 
structures.  These  problems  make  use  of  refined  rectangular  and  tri¬ 
angular  finite  elements.  Among  the  problems  considered  are  rectangu¬ 
lar,  circular,  and  annular  plates  with  various  boundary  conditions. 
Once  again,  comparisons  with  results  of  other  available  analyses 
are  favorable. 
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Problems  of  combined  bending  and  stretching  of  plates  are  also  con¬ 
sidered.  Results  are  obtained  for  rectangular  and  circular  plates. 
Results  for  combined  geometric  and  material  nonlinearity  are  pre¬ 
sented  for  beams  and  arches. 

55.  Mallett,  R.  H.,  AUTOMATED  METHOD  FOR  THE  LARGE  DEFLECTION  AND  IN¬ 
STABILITY  ANALYSIS  OF  THREE-DIMENSIONAL  TRUSS  AND  FRAME  ASSEMBLIES, 
AFFDL-TR-66-102,  December  1966. 

The  computer  program  presented  in  this  report  was  developed  to  pre¬ 
dict  large  deflection  behavior  of  three-dimensional  truss  and  frame 
assemblies.  The  solutions  are  obtained  by  the  direct  minimization  of 
the  total  potential  energy  with  respect  to  the  displacement  variables 
rather  than  by  solving  nonlinear  matrix  equations.  Sample  problems 
are  presented  to  demonstrate  the  analysis  capability  developed.  In¬ 
structions  for  the  preparation  of  the  input  data  and  the  FORTRAN  IV 
source  program  listing  are  included. 

56.  Przemieniecki,  J.S.,  MATRIX  METHODS  IN  STRUCTURAL  MECHANICS,  AFFDL- 
TR-66-80,  Conference  held  October  26-28,  1965. 

The  Conference  on  Matrix  Methods  in  Structural  Mechanics  held  at 
Wright-Patterson  Air  Force  Base  on  26-28  October  1965  was  sponsored 
jointly  by  the  Air  Force  Flight  Dynamics  Laboratory,  Research  and 
Technology  Division,  Air  Force  Systems  Command,  and  the  Air  Force 
Institute  of  Technology,  Air  University.  The  purpose  of  the  confer¬ 
ence  was  to  discuss  the  recent  developments  in  the  field  of  matrix 
methods  of  structural  analysis  and  design  of  aerospace  vehicles. 

57*  Berke,  L.,  PROCEEDINGS  OF  THE  SECOND  CONFERENCE  ON  MATRIX  METHODS  IN 
STRUCTURAL  MECHANICS,  AFFDL-TR-68-150,  Conference  held  October  15-17, 
1968. 

The  Second  Conference  on  Matrix  Methods  in  Structural  Mechanics,  spon¬ 
sored  by  the  Air  Force  Institute  of  Technology  (AFIT),  Air  University, 
and  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL),  Air  Force 
Systems  Command,  was  held  on  15-17  October  1968.  The  purpose  of  the 
conference  is  to  discuss  the  recent  developments  in  matrix  structural 
analysis  and  design  of  structural  systems.  This  volume  contains  all 
the  papers  presented  at  the  conference. 

Forty  papers  were  presented  at  the  conference  in  seven  sessions,  en¬ 
titled  Structural  Height  Optimization,  Dynamics,  General  Elements, 
Curved  Elements,  Applications,  General  Methods,  and  Nonlinear  Analy¬ 
sis.  The  papers  covered  practically  all  major  aspects  of  recent  re¬ 
search  and  development  work  on  matrix  methods  in  structural  mechanics. 

58.  Syraonds,  P.S.,  SURVEY  OF  TECHNICAL  METHODS  OF  ANALYSIS  FOR  PLASTIC 
DEFORMATION  OF  STRUCTURES  UNDER  DYNAMIC  LOADING,  BU/NSRDC/l-67, 

Brown  University,  Providence,  R.I.,  1967. 
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This  survey  attempts  to  make  a  critical  study  of  methods  described  in 
the  literature  for  the  analysis  of  metal  structures  under  dynamic 
loading  to  plastic  deformation. 

Analyses  have  now  appeared  in  the  literature  of  a  considerable  variety 
of  structures  of  engineering  interest.  They  include  beams  (under  many 
conditions  of  loading,  support,  and  materials),  rings,  arches ,  frames, 
(simple  rectangular  bents),  plates  (circular  and  rectangular),  mem¬ 
branes  (i.e.,  plates  with  deflections  greatly  exceeding  the  thickness), 
and  shells  (axially  symmetric  loading  on  cylinders,  spheres,  and 
spherical  caps).  Most  of  these  have  been  obtained  by  a  rigid-plastic 
type  of  analysis  (in  which  strain  rates  are  assumed  zero  unless  a 
yield  condition  is  satisfied).  A  few  have  been  obtained  by  wholly 
numerical  approaches  of  finite-difference  type. 

Experiments  reported  in  the  literature  have,  in  most  cases,  shown  that 
the  actual  permanent  deflections  are  smaller  than  those  predicted  on 
the  basis  of  plastic  properties  determined  by  quasi-static  tests,  the 
predictions  often  being  in  error  by  as  much  as  100  percent  or  more  for 
mild  steel,  with  smaller  discrepancies  for  other  metals,  such  as 
aluminum  alloys  or  high-strength  steels.  Strengthening  under  condi¬ 
tions  of  rapid  straining  has  been  considered  the  principal  cause  of 
such  discrepancies;  when  it  has  been  possible  to  modify  the  analysis 
to  take  account  of  the  increase  of  yield  and  flow  stresses  at  high 
strain  rates,  much  better  agreement  has,  in  most  cases, been  obtained. 

59«  Semonian,  J.W.,  and  Anderson,  P.A.,  AN  ANALYSIS  OF  THE  STABILITY  AND 
ULTIMATE  BENDING  STRENGTH  OF  MULTIWEB  BEAMS  WITH  FORMED  CHANNEL  WEBS, 
NACA  Technical  Note  TN  3232,  1954. 

Design  curves  and  procedures  are  presented  for  calculating  the  stress¬ 
es  for  instability  and  failure  of  multiveb  beams  with  formed-channel 
webs.  The  ultimate  bending  strength  of  this  type  of  construction  is 
shown  to  depend  upon  the  deflectional  stiffness  of  the  web  attachment 
flanges.  A  simple  criterion  is  also  given  for  predicting  whether  a 
multiweb  beam  with  a  given  attachment-flange  design  will  be  suscept¬ 
ible  to  a  wrinkling  instability  or  will  buckle  as  if  the  webs  are 
integrally  joined  to  the  cover  skins. 

The  criteria  for  predicting  buckling  and  failure  stresses  are  com¬ 
pared  with  experimental  data.  These  criteria  are  sensitive  to  the 
offset,  pitch,  and  diameter  of  the  rivets  used  on  the  web  attachment 
flanges,  and  the  riveting  specification  is,  therefore,  emphasized  as 
an  important  design  consideration. 

60.  Semonian,  J.W.,  and  Peterson,  J.P.,  AN  ANALYSIS  OF  THE  STABILITY  AND 
ULTIMATE  COMPRESSIVE  STRENGTH  OF  SHORT  SHEET -STRINGER  PANELS  WITH 
SPECIAL  REFERENCE  TO  THE  INFLUENCE  OF  THE  RIVETED  CONNECTION  BETWEEN 
SHEET  AND  STRINGER,  NACA  Technical  Report  TR  1255,  1956. 
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TEST  DATA 


DYNAMIC  TEST  INSTALLATION 

The  dynamic  test  setup  uses  existing  components.  The  drop  carriage  assem¬ 
bly  is  the  same  one  that  v?.s  utilized  during  the  drop  test  of  the  fuselage 
bumper,  the  results  of  which  are  reported  in  Reference  1.  Between  the 
test  specimen  support  and  the  ground  are  installed  six  load  cells  (two 
each  along  the  north  and  south  edges  of  the  support  and  one  each  along 
the  east  and  west  edge  of  the  support) .  The  load  cells  are  used  for  tests 
7  and  8.  In  the  previous  dynamic  tests  (4,  5,  6),  the  specimen  support 
is  grouted  to  the  concrete  slab  ground.  The  installation  provides  suffi¬ 
cient  free-fall  clearance  (l4  ft)  to  perform  30-ft/sec  impact  tests.  The 
setup  has  the  flexibility  of  performing  higher  impact  velocity  drop  tests 
by  adding  additional  frames,  thus  increasing  the  available  free-fall  dis¬ 
tance.  Figure  1  shows  a  layout  of  the  test  installation  and  notes  the 
various  items  and  assemblies  that  form  a  part  of  the  complete  installa¬ 
tion. 

TEST  SPECIMENS 

A  complete  description  of  the  test  specimen  overall  dimensions,  bulkhead 
and  end  beam  web  thicknesses,  number  and  thickness  of  angles,  design  con¬ 
figuration,  nominal  weight,  and  type  of  test  performed  for  each  specimen 
is  shown  in  Volume  I,  Table  XVII.  A  complete  set  of  drawings  for  the 
various  design  configurations  is  presented  herein  as  described  in  Table 
II  below: 


Table  II.  SPECIMEN  DRAWING  IDENTIFICATION 

Figure 

Applicable 

Specimen 

Drawing 

X-16628A 

2 

1 

-1 

3 

2,3,4 

-2 

4 

5,6 

-3 

5 

7,8,9 

-4 

6 

10,11,12 

-5 

RECORDED  SCAN'  (TIME)  HISTORIES 

A  complete  set  of  pertinent  recorded  test  data  is  presented  in  Figures  7 
through  80.  Load -deflection  curves  for  each  of  the  test  specimens  are 
presented  in  Volume  I  in  the  section  entitled  Substructure  Test  Program. 
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Figure  2.  Drawing  for  Specimen  1. 
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Figure  3-  Drawing  for  Specimens  2,  3>  and  4. 
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Figure  6.  Drawing  for  Specimens  10,  11,  and  12. 
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The  data  in  this  section  includes  load,  deflection,  acceleration  (tests 
4-8)  and  strain  versus  scan  plots.  Table  III  shovrs  the  test  data  available 
in  the  section  for  each- specimen. 

The  sampling  rate  for  test  4  is  750  scans  /second.  The  sampling  rate  for 
tests  5  through  8  is  1500  scans/second;  thus  for  the  dynamic  tests  the 
time  in  seconds  can  be  obtained  by  dividing  the  abscissa  scale  (scan)  by 
750  for  test  4  and  by  1500  for  tests  5  through  8.  Positive  strain  is  com¬ 
pression  for  tests  1,  2,  3>  9»  10,  H»  12  and  tension  for  tests  4,  5»  6>  7 
and  8. 
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Figure  15.  Deflection  Versus  Scan,  Test  Specimen  2. 
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Figure  17.  Strain  Gage  2  Versus  Scan,  Test  Specimen  2. 


Figure  25.  Strain  Gage  2  Versus  Scan  Test  Specimen  3. 
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Figure  27.  Strain  Gage  5  Versus  Scan,  Test  Specimen  4 


Figure  ^0,  Deflection  Versus  Scan,  Test  Specimen  5 
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Figure  59.  Acceleration  (Center)  Versus  Scan,  Test  Specimen  7 
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Figure  t>0,  Accel?ra‘ ion  (.'<e:-i'  Versus  Scan,  Test  Specimen 
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Figure  73*  Deflection  Versus  Scan,  Test  Specimen  10. 
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PROGRAM  "KRASH”  REFINEMENT 


GENERAL 

Program  KRASH  was  developed  as  an  analytical  tool  to  be  used  during  prelim¬ 
inary  phases  of  design.  In  Reference  1  the  results  of  a  parametric  study 
were  integrated  into  an  iterative  design  procedure  by  which  a  tradeoff 
between  potential  incremental  cost  and/or  weight  versus  incremental  improve¬ 
ments  in  crashworthiness  capability  could  be  accomplished.  The  results  of 
that  study  and  a  subsequent  study  for  the  U. S.  Army  (Reference  106)  showed 
the  potential  benefits  that  could  be  achieved  for  designers  using  a  com¬ 
prehensive  but  unsophisticated  analytical  approach.  However,  as  is  the 
situation  with  rapidly  developed  advances  in  the  state-of-the-art  analyti¬ 
cal  techniques,  there  arc  generally  areas  which  can  be  simplified  for  users 
not  thoroughly  familiar  with  the  approach. 


In  an  effort  to  facilitate  a  designer's  usage  of  program  KRASH  and  to  in¬ 
corporate  the  results  of  the  literature  survey,  load  sensitivity  study,  sub¬ 
structure  analysis,  and  tests,  the  program  was  revised.  In  particular,  the 
input  data  format  was  changed  for  ease  of  data  input  and  subsequent  data 
changes,  as  would  be  required  during  parameter  tradeoff  studies.  In  addi¬ 
tion,  flexibility  was  added  by  the  manner  in  which  the  Stiffness  Reduction 
factors  (KR's)  are  input.  The  revised  program  was  run  for  the  following 
condj  tions : 

•  Correlation  case  31-52  from  Reference  1  (to  demonstrate  that  the 
revised  program  format  was  compatible  with  previous  results:  23- 
ft/sec  vertical  velocity  combined  with  l8.5-ft/sec  lateral  velocity 
impact) . 

•  Three-dimensional  velocity  (~  40  ft/sec  combined  velocity  impact  - 
27.5  ft/sec  longitudinal,  23  ft/sec  vertical,  18.5  ft/sec  lateral). 

•  Upper  mass  penetration  into  a  specified  occupiable  volume. 

•  Simplified  blade  contact. 

•  Utilization  of  load-deflection  data  obtained  from  the  12  substruc¬ 
ture  tests,  performed  during  the  study,  and  related  to  actual  fuse¬ 
lage  structure  size. 

PROGRAM  " KRASH"  INPUT  FORMAT  REVISIONS 

The  input  format  changes  are  divided  into  three  categories: 

•  Reordered  data. 

•  Standardization  of  certain  inputs. 

•  Allowance  for  more  general  KR  curves. 
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First,  the  input  data  is  rearranged  so  that  all  mass  associated  data 
(subscript  i)  is  in  one  block,  followed  by  all  external  spring  data  (sub¬ 
script  ik)  and  then  all  internal  beam  data  (subscript  ij).  Next,  standard 
values  are  assigned  to  certain  input  data  items  unless  otherwise  specified. 
Thus,  much  repetitive  data  input  is  eliminated .  The  following  quantities 
are  given  standard  values  unless  otherwise  specified: 


Quantity 

. 

Symbol 

Standard 

Value 

Angular  Momenta 

Hei 

0 

Euler  Angles 

±tt  Q  ••  A" 

4>  ±,  ▼  i 

0 

Aerodynamic  Lift 

lc. 

1 

0 

Stiffness  Reduction  Factors 

KR.  . 
i.1 

1 

Failure  Deflection 

Uraaxiji 

100 

Thus,  the  angular  momenta  of  the  masses  (He^ )  and  the  aerodynamic  lift  (lci) 

are  normally  zero,  as  are  the  Euler  angles  6",  9”,  relating  mass  body- 

fixed  axes  to  airplane  c.g.  axes.  For  linear  internal  beam  elements,  KR  =1 
for  the  entire  run,  so  only  nonlinear  KR  data  need  be  input.  Similarly, 
v  „  ,  need  only  be  input  for  those  elements  where  rupture  is  expected  to 

occur;  for  elements  not  specified,  a  failure  deflection  of  100  inches  and 
rotation  of  100  radians  are  assumed.  Thus,  the  nonspecified  beams  will  not 
rupture  during  the  run. 

The  input  has  also  been  revised  to  allow  more  general  KR  (stiffness  reduc¬ 
tion  factor)  curves.  Previously,  each  KR  curve  could  have  only  six  data 
points,  equally  spaced.  Now  each  curve  can  have  up  to  15  data  points,  with 
any  desired  sr,-~ing.  This  greatly  facilitates  modeling  more  complex 
load-deflection  relationships.  However,  the  load-sensitivity  study  results 
indicate  that  for  many  structural  elements  a  very  simple  KR  representation 
is  adequate.  The  KR  representation  for  the  various  load-deflection 
categories  is  discussed  in  Volume  I  (Design  Procedures). 

The  relationship  between  the  number  of  masses,  internal  beam  elements-  and 
KR  tables  is  governed  by  the  following  expression: 

Total  number  of  bytes  =  235,136  +  1244  x  number  of  lumped  masses  (n) 

+  1372  x  Internal  Beam  Elements  (IBE)  +  19B  x  number  of  KR  tables 

e.g.  N  =  80 
IBE  =  100 
KR  =  120 

bytes  =  235136  +  260, 480  =  495616  =  484k 
where  IK  =  1024  bytes. 
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Thus ,  program  KRASH  has  sane  flexibility  with  regard  to  the  size  of  problems 
to  which  it  can  be  applied.  Depending  on  the  type  of  vehicle  that  is  being 
analyzed  the  program  can  be  redimensionali zed  to  provide  the  most  practical 
and  economical  solution. 

The  three-dimensional  impact  mass  penetration,  sinplified  rotor  blade  con¬ 
tact  and  updated  load-deflection  data  cases  are  described  in  Volm.'e  I  under 
the  section  entitled  Program  KRASH  Refinement. 

ENERGY  BALANCE  EQUATIONS 

The  primary  objective  ox  a  crash  analysis  in  the  preliminary  design  phase 
is  to  determine  how  to  absorb  the  initial  vehicle  kinetic  energy  while 
maintaining  a  livable  environment  for  the  occupants.  This  task  is  greatly 
facilitated  if  information  regarding  the  spatial  distribution  of  the  energy 
flow  through  the  vehicle  is  available.  With  the  objective  in  mind,  energy 
balance  equations  are  developed.  These  equations  do  not  alter  the  existing 
computational  procedures  or  results;  they  merely  provide  additional  infor¬ 
mation  to  assist  in  understanding  how  the  initial  kinetic  energy  is 
absorbed. 

The  total  system  energy  at  any  time  is  given  by  the  following  expression: 

Etqt  =  KE  +  PE  +  SE  +  DE  +  CE  (l) 

w*iere  Etqt  =  Total  system  energy 

KE  =  Total  kinetic  energy 

PE  =  Total  potential  energy 

SE  =  Total  strain  energy  absorbed 

DE  =  Total  damping  energy  dissipated 

CE  =  Total  crash  spring  (external  spring)  energy 
absorbed . 

The  total  system  energy  E^  remains  constant  during  the  analysis.  The 

total  kinetic  and  potential  energies  for  each  mass  over  the  number  of 
masses  (N)  are  shown  in  Equation  (2). 

N  N 

KE  =  2  KE.  ?E  =  E 

i  =  1  i  =  1 

The  total  strain  and  damping  energies  are  obtained  by 
and  danping  energy  for  each  interval  beam  element  (ij 
pairs . 


(2) 


summing  the  strain 
pair)  over  the  M  ij 
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DE  = 
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DE 
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(3) 


SE^  results  from  the  elastic-plastic  behavior  of  the  beam,  and  DE^ 

results  from  its  damping  properties.  The  total  crash  spring  energy  re¬ 
sults  from  summing  the  energies  for  all  the  individual  crash  springs  (ik 
pairs)  over  all  the  P  ik  pairs. 


P 

CE  =  £ 


00 


ik  =  1 


CE., 
ik 


Referring  to  Equation  (l),  at  time  zero  all  energies  are  zero  except  KE, 

PE  and  The  potential  energy  is  referenced  to  the  ground  plane. 

After  impact  with  the  ground,  the  kinetic  and  potential  energies  decrease; 
damping, strain  and  crash  spring  energies  all  increase  to  keep  ETOTconstant. 

In  Equation  (3),  the  summations  over  the  ij  pairs  exclude  those  ij  pairs 
that  are  identified  in  the  input  as  DRI  elements,  which  are  described  in 
Reference  1,  page  73 •  This  is  done  because  these  ij  beams  and  their 
masses  are  isolated  from  the  rest  of  the  system;  the  forces  in  DRI  beam 
ij  drive  mass  j  but  not  mass  i.  Also,  the  summations  in  Equation  (2) 
for  the  kinetic  and  potential  energies  exclude  mass  j  in  a  DRI  ij  pair, 
since  this  mass  is  isolated  from  the  system.  Thus,  if  ij  pairs  6-9  and 
11-15  are  defined  as  DRI  elements,  the  summations  for  SE  and  DE  will 
exclude  these  ij  pairs  and  the  summations  for  KE  and  PS  will  exclude 
masses  9  and  15. 

The  kinetic  energy  for  each  mass,  including  translational  and  rotational 
components,  is  simply 
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where 


and 
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1  q. 
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|Mi|  is  the  6x6  inertia  matrix  for  mass  i,  and  jvel^j  a  six-element 

vector  of  the  linear  and  angular  velocity  components  of  mass  i  in  body- 
fixed  axes. 

The  i  mass  potential  energy,  referenced  to  the  ground  plane,  is  given  by 


PEi  =  wi  zi 


(6) 


Note  that  z ^  is  positive  downward,  measured  from  the  ground  plane. 


The  strain  energy  in  internal  beam  ij  is  computed  as  a  continuous  summa¬ 
tion  of  the  incremental  energy  contributions  from  each  integration  interval. 
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F. .  Avb  is  the  internal  beam  force  (or  moment)  in  the  1 


th 


(7) 


direction, 
,th 


multiplied  by  the  incremental  beam  deflection  (or  rotation)  in  the  i 
direction;  this  is  the  incremental  strain  energy  for  the  integration  inter¬ 
val  being  considered.  This  straightforward  formulation  automatically  ac¬ 
counts  for  the  complexities  of  nonlinear,  coupled  deflections  and  unloading - 
reloading  behavior,  since  these  are  considered  in  the  calculation  of  F. .  . 


The  damping  energy  dissipated  in  internal  beam  ij  is  computed  in  a  Similar 
manner; 


M= 


urrent 
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FD. ..  is  the  internal  beam  damping  force  (or  moment)  in  the  £ 
ijf 


th 


(8) 

direction. 
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The  crash  spring  energy  is  also  computed  as  a  summation  of  incremental 
energy  changes.  The  crash  spring  energy  resulti:ig  from  all  the  ik  springs 
attached  to  mass  i  is 


current 


previous 


XCiYCi-’  NCi 


(9) 


X^,  YCi’  *'■’  NCi  are  S^JC  forces  ant*  momen'ts  acting  on  mass  i  in 

body-fixed  axes,  resulting  from  all  ik  external  springs  attached  to 
mass  i.  These  are  given  by  Equations  (66)  in  Reference  1,  page  6l.  The 


vector  ax 'p  Ay^,  Az'^,  Ainp.^,  Ainr^J  is  made  up  of  the  six  incremental 

deflections  and  rotations  of  mass  m  in  the  same  body-fixed  axes.  The 
first  three  terms  of  this  vector  aregiven  by  a  simple  rotation  transfor¬ 
mation  of  the  incremental  deflections  in  ground  axes,  which  are  obtained 
directly  from  numerical  integration  of  the  equations  of  motion. 


(10) 


The  last  three  terms  of  the  incremental  displacement  -vector  in  Equation  (9) 
are  the  incremental  changes  in  the  integrals  of  the  angular  velocities  of 
mass  un,  p^,  q^,  r^.  These  are  the  incremental  rotations  of  mass  m^  in 

body -fixed  axes  for  the  integration  interval  being  considered. 

The  negative  sign  in  Equation  (9)  results  from  the  fact  that  the  forces 
acting  on  mass  m^  rather  than  the  forces  within  the  ik  springs  sure 
being  considered.  A  positive  deflection  of  spring  ik  results  in  a 
negative  force  on  mass  m. .  The  energy  calculated  in  Equation  (9) 

includes  tne  energy  dissipated  by  the  sliding  of  spring  ik  on  the 
ground  with  a  friction  coefficient ;  hence  the  forces  in  that  equation 
include  the  ground  drag  loads  due  to  friction. 

The  crash  spring  energies  in  Equation  (9)  are  not  yet  in  the  form  desired; 
the  crash  spring  energy  due  to  each  spring  ik  must  be  separated  out  of  the 
total  crash  spring  energy  associated  with  mass  nu  ,  CE^.  This  is  done 
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simply  by  substituting  Equations  (66)  of  Reference  1  for  the 

jxci,  Yci,  Ncij  terms  in  Equation  (9).  This  reformulates  Equation 

(9)  into  a  function  of  the  individual  crash  spring  forces  FSP^,  where  i 

and  k  refer  to  the  ik  spring  and  j  refers  to  the  direction  of  the 
forces  on  the  ik  spring.  These  forces  are  shown  in  Figure  10  on  page  62 
of  Reference  1.  The  final  equation  for  the  crash  spring  energy  associated 
with  each  spring  ik  is  the  following: 


current 


previous 


3 

S  FSPi3J[  Avc^  +  KEMik  (11) 
J  1 


where 


and 


TEBMik  -  FSP131  JU  41.^  -  FSPia  ^  Ain^  (k  =  1) 

™»ik  =-FSP132  Ainp,  *  FSPil2  |.2  Ainr.  (k  -  2) 

raRMik  *  FSPi23  *13  MnPi  ’  FSPil3  *13  *lmli  (k  ■  3) 


This  is  the  crash  spring  energy  to  be  used  in  Equation  (4)  to  obtain  the 
total  crash  spring  energy,  CE,  for  the  entire  vehicle.  Note  that  CE  could 
have  been  determined  more  directly  by  summing  the  CE^  from  Equation  (9) 
over  all  the  masses  nu  .  However,  the  CE^  associated  with  each  spring  ik 

is  calculated  so  that  it  can  be  determined  how  much  each  external  spring 
is  contributing  to  CE. 

ENERGY  BALANCE  DATA 


The  energy  balance  technique  described  in  the  preceding  section  and  the 
information  presented  in  the  following  paragraphs  describe  a  type  of  format 
which  can  be  of  benefit  to  a  designer  during  the  preliminary  stage  of  de¬ 
sign. 
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The  output  format  for  the  energy  balance  data  is  shovm  in  Figures  8l,  82 
and  83.  For  a  condition  of  42 -ft/sec  vertical  velocity  impact  and  using 
program  KRASH  with  the  existing  UH-1H  helicopter  math  model  consisting  of 
31  masses  and  38  internal  beam  elements,  Figure  8l  shows  that  at  time  * 

0.0  second  only  contributions  to  the  energy  term  are  kinetic  and  potential 
energies.  The  kinetic  energy  represent,’  the  l/2  M\r  terms  while  the  po¬ 
tential  energy  is  a  function  of  the  height  the  respective  masses  are  above 
the  ground  at  the  initiation  of  impact.  Hie  internal  beam  energy  contri¬ 
butions  due  to  strain  and  damping  are  zero  at  this  time,  as  is  the  exter¬ 
nal  spring  crushing  energy .  Figure  82  shows  the  energy  data  output  during 
the  run  at  time  =  .042  second  after  impact.  At  this  time  kinetic  energy 
has  decreased  from  85.85$  of  the  total  energy  to  58.6$  of  the  total  energy. 
Similarly,  the  potential  energy  has  decreased  from  14.15$  to  8.96$  of  the 
total  energy.  At  this  time  crushing  energy,  strain  energy  and  damping 
energy  as  a  percent  of  the  total  energy  have  increased  to  9.8$,  7.5$  and 
15.2$,  respectively.  Furthermore,  Figure  82  snows  the  following: 

•  distribution  of  kinetic  energy  and  potential  energy  by  mass  item 

•  distribution  of  strain  energy  and  damping  energy  by  beam  element 

•  distribution  of  crushing  energy  by  external  spring  element 

The  negative  terms  in  the  damping  energy  output  for  beam  elements  2  and  4 
are  of  the  order  of  less  than  .0003  percent  of  the  total  energy.  The 
values  are  obviously  insignificant  and  are  considered  to  reflect  damping 
inherent  in  the  integration  routine. 

Figure  83  shows  a  summary  of  the  energy  data  which  is  presented  at  the  end 
of  the  computer  run.  The  summary  shows  for  each  time  increment  the  total 
energy  and  the  amount  and  percentage  contribution  from  kinetic,  potential, 
strain  damping  and  crushing.  The  data  in  the  summary  shown  in  Figure  83  * 
indicates  that  for  the  illustrated  case  the  vehicle  rebounds  at  approxi- 
mateiy  .064  second  after  impact.  The  kinetic  energy  reduces  from  impact 
until  .064  second  after  impact, then  starts  to  increase  again.  The  summary 
data  presented  in  Figure  83  also  shows  that  crushing  energy  reaches  a  max¬ 
imum  of  51.73$  of  the  total  energy  at  ,062  second  after  impact  and  that 
strain  en.rgy  reaches  a  maximum  of  21.93$  of  the  total  energy  at  .076 
second  aftsr  impact. 
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Figure  8l.  Energy  Balance  Data  Output  at  t  =  0.0  Second 
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Figure  82.  Energy  Balance  Data  Output  at  t  =  0.0^2  Secooi. 
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This  section  describes  all  the  input  necessary  to  run  the  program,  and  the 
output  available  from  the  program. 


Input 


The  input  data  format  is  described  in  detail  in  this  section  and  is.  shown 
in  Figure  84.  Unless  otherwise  specified,  all  quantities  are  input  in 
inch,  pound,  second,  and  radian  units.  For  cards  0102  and  following  the 
input  format  is  6E12.0  unless  otherwise  noted.  Each  card  has  6  fields; 
each  field  is  12  columns  wide.  As  an  example  of  this  format,  the  number 
126.08  can  be  input  in  the  following  ways: 


Blank  columns  are  treated  as  zeros.  When  the  E  format  is  used,  the  expo¬ 
nent  must  be  right  justified  in  the  field.  Integer  formats  with  field 
widths  of  1,  2  or  3  are  also  used.  Ihe  numbers  in  these  locations  must  be 
right  justified  integers.  Sequence  numbers  in  columns  77  through  80  should 
be  used  corresponding  to  those  shown  in  the  input  format  to  facilitate 
deck  assembly  and  changes.  All  tabular  input  is  linearly  interpolated  be¬ 
tween  input  values  and  extrapolated  beyond  the  two  end  values,  if  neces¬ 
sary. 


Card  0100  -  This  card  contains  the  title  for  the  case  being  analyzed. 
All  text  material  on  card  0100  is  reproduced  at  the  top  of  every  page 
of  the  output  and  on  every  plot. 


Card  0101  -  N  is  the  total  number  of  lumped  masses.  The  maximum  allow¬ 
able  number  of  masses  is  80.  APrint/At  is  the  integer  multiple  of  At 
at  which  output  is  printed.  At  is  the  numerical  integration  time  inter- 
val.  is  the  time  span  being  analyzed. 

-  These  cards  contain  the  overall  vehicle  ini¬ 
tial  conditions,  x^,  yQ,  and  zQ  are  the  ground  axes  components  of  the 

initial  c.g.  velocity,  p',  q',  and  r'  are  the  c.g.  coordinate  system 
components  of  the  initial  angular  velocity  of  the  vehicle,  p*  is  the 
roll  velocity,  q'  the  pitch  velocity,  and  r'  the  yaw  velocity.  9',  0', 
and  i|f 1  are  the  Euler  angles  relating  the  initial  position  of  the  vehi¬ 
cle  to  ground  coordinates.  4>;  is  the  roll  angle,  0'  the  pitch  angle, 
and  f  the  yaw  angle.  zQ  is  the  negative  of  the  initial  vehicle  c.g. 
height  above  ground.  If  this  input  is  zero  (blank),  the  initial  condi¬ 
tion  subroutine  computes  a  value  of  z~  so  that  the  lowest  extremity  of 
the ’vehicle  is  .1  inch  above  the  ground. 


Cards  0102 
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Ca  ds  0201  through  02 XX  -  N  of  these  cards  are  used  to  input  the  weights 
{ncc  masses)  of  the  N  lumped  masses. 


Cards  0301  through  03XX  -  N  of  these  cards  sure  required.  Each  card 
inputs  the  six  moments  and  products  of  inertia  for  the  ith  maxx,  i  =  1, 

<?,  .  •  • ,  N. 

Cards  C401  through  04XX  -  N  of  these  cards  are  used  to  input  the  coordi¬ 
nates  of  the  N  lumped  masses,  xj  is  the  Fuselage  Station  (increasing 
aft),  y"  is  the  Butt  Line  (positive  left),  and  zV  is  the  Water  Line 
(increasing  upward). 

Card  0500  -  NI  is  the  number  of  masses  having  nonzero  Hexi,  He  ,  Hezi 

(angular  momenta)  or  4>",  0^,  which  are  the  Euler  angles  relating  the 

c.g.  axes  to  the  ith  mass  body  fixed  axes.  The  above  quantities  are 
normally  zero,  so  only  nonzero  values  are  input  on  cards  0501  through 
05NI.  i.,  i2,  ...,  iKI  are  the  actual  i's  or  mass  numbers  that  have  non¬ 

zero  input  data  for  any  of  the  above  quantities.  NI  and  i  through  i^ 

are  input  in  integer  format  12.  If  NI  equals  zero,  cards  0501  through 

05NI  are  not  input;  however,  card  0500  is  always  required.  If  NI  equals 

zero,  a  bland  card  for  card  0500  is  input. 

Cards  0501  through  05NI  -  As  noted  above,  these  cards  are  used  to  inpu: 
any  nonzero  values  of  He^,  He^  or  6^,  0l\  If”,  One  card  is  used  for 

each  mass  i  having  nonzero  input.  The  masses  are  ordered  according  to 
the  sequence  specified  on  card  0500.  He^^,  He^,  and  He^  are  the  body 

axes  components  of  the  angular  momenta  of  masses  m^,  d.;e  to  rotation  of 
internal  masses  within  n^.  These  are  normally  zero.  oV,  8”,  are 

the  Euler  angles  from  the  c.g.  axes  to  the  ith  mass  axes.  If  the  ith 
mass  body-fixed  axes  are  parallel  to  the  vehicle  c.g.  coordinate  axe3. 
which  is  usually  the  case,  these  are  all  zero.  Note  that  if  any  non¬ 
zero  values  are  input,  then  8.  .  and  <|r.  .  on  cards  0901  through  090M  must 
be  input.  1J 

Cards  0600  through  06nI  -  NI  is  the  number  of  masses  having  nonzero  lCi, 

which  are  the  aerodynamic  lift  constants  used  in  Equation  (5)  on  page  30 
of  Reference  1.  ij.  ig»  ...»  i^j  sore  the  actual  mass  numbers  having  non¬ 
zero  lc^:  If  all  1^  are  zero,  card  0600  is  input  as  a  blank  card  and 

cards  06d  through  06NI  are  not  input.  One  card  is  used  for  each  non¬ 
zero  lCi  input,  with  the  ordering  defined  on  card  0600. 

Cards  0701  through  070p  -  As  many  cards  as  necessary  are  used  to  input 
data  for  the  external  springs.  The  maximum  allowable  number  of  external 
springs  is  50.  Each  spring  is  identified  by  an  ik  pair,  where  i  defines 
the  mass  to  which  the  spring  is  attached  and  k  defines  the  direction  of 
the  spring,  k  =  1,  2  and  3  correspond_to  springs  in  the  x,  y  and  z  di¬ 
rections,  in  mass  i  body-fixed  axes.  is  the  free  length  of  the  ex¬ 
ternal  spring  ik.  is  positive  if  it  radiates  out  from  mass  nu  m 


■  sV  „ 


the  positive  direction  of  the  ith  mass  body  axes;  1^  is  negative  if  it 

radiates  in  the  opposite  direction.  Springs  in  both  directions  are  not 
allowed.  mu^  is  the  friction  coefficient  between  the  ground  and  the 

end  of  the  ik  spring,  is  the  linear  unloading  stiffness  and  also 

the  bottoming  stiffness  for  the  ik  spring. 

If  no  external  springs  are  used,  card  0701  is  input  blank  and  the  next 
card  input  is  0901. 

Card  07XX  -  A  blank  card  in  this  location  is  required  to  terminate  the 
above  ik  external  spring  list. 

Cards  0801  through  080p  -  External  spring  load-deflection  curve  param¬ 
eters  are  input  on  t^.ese  cards,  one  card  per  ik  spring,  ordered  as  in 
the  ik  list  on  cards  0701  through  070p.  The  program  is  written  so  that 
a  table  of  the  following  form  is  input:  ' 


This  table  is  defined  by  six  parameters  Sm,  S^,  S Bik,  S^,  FSP0Iik, 
and  FSPQ^k,  which  are  input  on  cards  0801  through  080p. 


Cards  0901  through  090M  -  These  cards  contain  the  Euler  angles  4^,  0^, 
and  i»..  for  all  internal  beams  lj .  The  beam  interconnections  are'  defined 

by  the  i's  and  j's  input,  i  must  be  less  than  j,  but  there  Is  no  other 
requirement  on  the  oi dering  of  the  ij  pairs.  is  always  input;  0^ 

and  t| need  not  be  input  if  <t>£,  8'.'  and  yV  on  cards  0501  through  05NI  are 

all  zero.  In  the  latter  case,  8 .  .  and  <|r.  .  are  computed  in  initial  condi- 

tior.s,  <t>.  .  will  normally  be  zero.  The  maximum  allowable  number  of  in- 

ternal  beams  is  100. 
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Card  09XX  -  A  blank  card  is  required  here  to  terminate  the  ij  pair 
internal  beam  list. 

Cards  1001  through  1006  -  These  six  cards  are  used  to  input  the  6x6 
linear  stiffness  matrix  jx^j  for  the  first  ij  beam,  listed  on  card 

0901 Each  card  inputs  one  row  of  the  matrix. 

Cards  1007  through  1XXX  -  As  many  cards  as  necessary  are  used  to  input 
all  the  remaining  6  x  o  jK^j  matrices.  These  matrices  must  be  ordered 

the  same  as  the  ij  pairs  are  ordered  on  cards  0901  through  090M. 

Cards  2001  through  200M  -  These  cards  are  used  to  input  the  C^  for  all 
the  ij  pairs  defined  on  cards  0901  through  090M.  is  the  damping 

ratio  (damping/critical  damping)  for  the  isolated  system  consisting  of 

masses  m.  and  m  connected  by  beam  ij.  Values  of  C.  between  .01  and 
^  J  ^*3 

.05  are  generally  representative  of  the  structural  damping. 


Cards  2101  through  2l0q  -  These  cards  are  used  to  specify  which  beam 
elements  (ij  pairs)  and  directions  (l)  have  nonlinear  stiffness  prop¬ 
erties.  1  varies  from  1  to  6  corresponding  to  the  x,  y,  z,  6,  9,  ijt, 
directions  in  beam  axes.  Only  beam  elements  and  directions  having  non¬ 
linear  properties  (KR  ^  l)  are  input.  For  each  ijl  combination  listed, 
a  table  of.  KR  (stiffness  reduction  factor)  versus  deflection  is  input 
on  the  following  cards.  The  NP's  are  the  number  of  points  in  the  fol¬ 
lowing  tables.  For  all  itjl  combinations  not  listed,  KR  =  1  for  the 
entire  run.  The  maximum  allowable  number  of  ijl  combinations  (non¬ 
linear  KR  tables)  is  80.  The  maximum  allowable  number  of  points  per 
table  is  15.  For  a  completely  linear  analysis,  card  2101  is  blank  and 
the  next  card  input  is  3001. 


Card  2 LXX  -  A  blank  card  is  required  here  to  terminate  the  ijl  list  of 
nonlinear  beam/direction  combinations. 

Cardo  2201  through  22NP^  -  These  NP^  cards  are  used  to  input  the  deflec¬ 
tion  (vb^)  versus  KR^  table  for  the  ijl  combination  listed  on  card 
2101.  KR. normally  starts  at  1.0  zero  deflection  (or  rotation),  cor- 

responding  to  linear  behavior,  and  varies  in  any  desired  manner  there¬ 
after  to  define  a  general  load-deflection  curve.  Note  that  KR  is  not 
the  actual  load,  but  rather  the  derivative  of  load  with  respect  to  de¬ 
flection,  i.e.,  the  slope  of  the  load-deflection  curve. 


Cards  22X1  through  2 XXX  -  As  many  cards  as  necessary  are  used  to  input 
the  KRj^  vs.  vb^  tables  for  the  remaining  ijl  combinations  listed  on 

cards  2100  through  210q  above.  Note  that  each  table  can  have  a  different 
number  .*f  data  points  NP. . 
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Cards  3001  through  300r  -  These  cards  sure  used  to  identify  which  internal 
beam  elements  ij  have  nonstandard  maximum  deflections  for  rupture 
(Vmaxijl) *  ®ie  3tandard  deflections  built  into  the  program  are  100 

inches  for  deflections  and  100  radians  for  rotations,  ihese  numbers 

were  deliberately  chosen  to  be  very  high  so  that  rupture  would  not  occur 

unless  reasonable  v  ,  were  specified.  If  no  nonstendard  v  „  ,  are 
max  s  max  s 

to  be  input,  card  3001  is  input  blank,  and  the  next  card  input  is  5001. 

Card  30XX  -  This  blank  card  is  required  to  terminate  the  preceding  ij 

list  for  nonstandard  v  .  . 

maxijl 

Cards  3101  through  310r  -  These  cards  are  used  to  input  the  nonstandard 
vmaxijl 's*  one  car<*  ^or  eac^  beam  specified  on  cards  3001  through 

300r.  Each  card  inputs  the  six  vmaxis  one  ij  beam,  ordered  x,  y,  z, 

<t>,  9,  t  in  beam  axes.  vmaxij]_  defines  the  maximum  allowable  deflection 

before  rupture  occurs.  After  rupture,  all  forces  in  beam  ij  go  to  zero, 
regardless  of  the  direction  in  which  the  rupture  occurred. 

Cards  5001  through  5025  -  These  cards  are  used  to  specify  the  time  his¬ 
tory  output  plots  desired.  The  only  input  for  these  cards  is  either  a 
1  or  a  blank.  A  1  results  in  the  output  of  a  time  history  plot  for  the 
response  quantity  indicated;  a  blank  results  in  no  plot  for  that  item. 

For  example,  a  1  in  the  13th  column  oi'  card  5003  specifies  that  a  time 
history  plot  of  z^  is  to  be  generated.  For  cards  5001  through  5012, 

the  column  number  in  which  the  1  is  input  indicates  which  mass  i  is  de¬ 
sired.  For  cards  5013  through  5024,  the  column  number  in  which  the  1  is 
input  indicates  which  ij  pair  is  desired,  where  the  ij  pairs  are  ordered 
as  on  cards  0200  through  02XX. 

Plots  are  available  for  the  displacements,  velocities  and  accelerations 
of  all  the  lumped  masses,  all  the  external  spring  compressions  S.^,  and 

all  the  beam  ij  total  deflections/rotations  (vb. .)  and  forces/moments 
(F,  .).  The  latter  two  items  are  in  beam  ij  axes?  Also  available  are 

plots  of  the  DRI  (Dynamic  Response  Index).  These  are  identified  by  the 
DRI  ij  element,  as  described  on  card  7000  below.  The  plot  variables  are 
labelled  automatically,  and  the  plot  scales  are  chosen  automatically. 

The  user  merely  has  to  specify'  which  plots  are  desired.  Up  to  150  plots 
can  be  requested  per  run.  Thirty  thousand  data  points  are  stored  for 
plotting,  with  the  plot  time  interval  depending  on  run  time  and  number 
of  plots  requested.  Thus,  if  the  maximum  of  150  plots  are  requested, 

200  points  will  be  saved  to  generate  each  plot.  If  the  maximum  run  time 
is  .2  second,  this  will  give  one  data  point  every  .001  second.  For  the 
types  of  problems  analyzed,  this  appears  to  be  a  marginally  acceptable 
resolution.  Requesting  fewer  plots  automatically  increases  bhe  number 
of  points  per  plot,  and,  hence,  the  resolution. 
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Card  6000  -  P  is  the  mass  to  be  used  to  locate  the  mass  penetration 
control  volume.  Hie  format  is  12.  If  no  mass  penetration  calculations 
are  required,  this  card  is  input  blank  and  card  6001  Ls  not  required. 

Card  6001  -  xn,  xp,  yn,  yp,  zn,  and  zp  are  the  six  distances  (all  posi¬ 
tive),  measured  from  the  control  mass  P  to  the  six  sides  of  the  control 
volume.  These  are  measured  along  the  positive  (p)  and  negative  (n) 
body-fixed  axes  for  mass  mp. 

Card  7000  -  This  card  is  used  to  specify  which  interned  beam  ij  elements, 
if  any,  are  to  be  treated  as  DRI  (Dynamic  Response  Index)  elements.  For 
example,  if  the  15th  ij  pair  listed  on  cards  0901  through  090M  is  to  be 
a  DRI  element,  then  a  1  is  input  in  column  15  of  card  7000.  For  DRI 
beam  elements,  the  forces  in  beam  ij  drive  mass  j  but  not  mass  i.  To 
understand  the  physical  significance  of  this,  refer  to  Reference  1,  page 
73,  which  explains  the  concept  of  a  DRI. 

Cards  8000  and  8001  -  Cards  8000  and  8001  are  used  to  specify  il  (mass- 
direction)  pairs  for  which  it  is  desired  to  input  a  time  history  table 
of  acceleration.  Thus,  any  mass  can  be  driven  by  a  specified  accelera¬ 
tion  in  any  direction.  The  6  directions  are  ordered  x,  y,  z,  6,  0,  f 
in  mass  i  body-fixed  axes.  The  format  is  24  (1211),  where  the  1  is 
right  justified  in  the  field  of  3  and  the  i  is  right  justified  in  the 
left  two  columns  of  the  field  of  3*  A  maximum  of  48  input  tables  are 
allowed.  For  example,  driving  mass  3  in  the  y  direction  would  be  speci¬ 
fied  as  |C32 | . 

When  any  acceleration  table  is  input,  that  acceleration  replaces  the 
normally  computed  acceleration,  so  it  becomes  a  driving  force  to  the 
system.  If  less  than  24  il  pairs  are  specified,  only  card  8000  is  re¬ 
quired  (omit  8001).  If  24  or  more  il  pairs  axe  input,  cards  8000  and 
8001  are  both  required  (if  exactly  24  il  pairs  are  input,  card  8001  will 
be  blank,  but  must  be  input). 

Card  8010  -  is  the  number  of  data  points  in  the  following  accelera¬ 
tion  time  history  table,  applicable  to  the  first  il  pair  on  card  8000. 

The  maximum  allowable  value  of  is  50. 

Cards  8011  through  801N^  -  cards  are  used  to  input  the  time  history 

table  of  acceleration  for  the  first  il  pair.  Each  card  inputs  one  time 
and  the  corresponding  acceleration. 


Cards  8XX0  through  8XXN^  -  The  remaining  acceleration  time  history 

tables  are  input  on  these  cards,  in  the  format  described  above.  The 
tables  are  ordered  according  to  the  sequence  of  il  pairs  defined  on 
cards  8000  and  8001. 
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First,  all  the  input  data  is  printed  out,  with  self-explanatory  identifying 
titles.  Next,  at  each  print  time  (=  APrint/At  x  At),  the  data  shown  in 
Figure  85  is  printed  out/  At  the  top  of  each  page,  the  case  identification 
data  is  printed  out  (from  input  card  0100 ).  Then  the  current  value  of  time 
t  is  printed. 

Next,  for  each  of  the  N  masses,  five  lines  of  data  are  output,  x.^,  y^,  and 

z^  are  the  ground  coordinates  of  m^.  x^,  y^,  and  zu  are  the  ground  axes 

components  of  the  velocity  of  m^.  u^,  v^,  are  the  time  derivatives  of 

Ui,  vi’  cnd  *  (^ese  61,6  not  equal  to  the  acceleration  of  nu,  as  can  be 

seen  from  Equation  (8l),  Reference  1.)  XACCEL,  YACCEL,  ZACCEL  are  the 
body  axes  components  of  the  acceleration  of  nn  . 

4>,  8,  and  ijr  are  the  Euler  angles  defining  the  attitude  cf  body  nn ,  in  roll, 
pitch  and  yaw.  4/ ,  8!!,  and  <|u  are  the  time  derivatives  of  <K ,  8^,  and  i|f. 
Pi,  q^j  and  ru  are  the  body  axes  components  of  the  angular  velocity  of  nn  ; 
they  are  the  velocities  in  roll,  pitch,  and  yaw,  respectively,  p^,  q^,  and 
r^  are  the  body  axes  components  of  the  tngular  acceleration  of  m^. 

Following  that  output,  the  running  time  sums  of  the  internal  forces  j* 

(Equation  (27),  Reference  l)  are  printed  out;  the  six  forces  and  moments 
for  each  ij  pair  are  printed  on  a  line  preceded  by  the  identifying  i  and  j. 
Next,  the  running  time  sums  of  the  beam  deflections  /y^ .  A  (Equation  (26), 
Reference  l)  are  printed  out  in  the  same  format.  I 

Finally,  the  external  spring  compressions  s^  (Equation  (47),  Reference  l) 

Eire  printed  out.  Each  line  starts  with  the  i,  followed  by  the 
k  =  l,  2,  3.  Only  values  of  i  for  which  a  spring  is  input  are  shown. 

During  the  course  of  the  run,  if  any  ruptures  or  control  volume  mass  pene¬ 
trations  occur,  the  appropriate  information  is  printed  out.  When  a  rupture 
occurs,  the  word  RUPTURE  is  printed,  followed  by  four  items: 

1.  The  ij  pair  that  ruptured,  where  the  numbering  corresponds  to  the 
ordering  of  the  ij  pairs  as  input  on  cards  0200  through  02XX. 

2„  The  1  (from  1  to  6)  indicating  in  which  direction  (in  beam  axes) 
the  rupture  occurred.  These  are  ordered  x,  y,  z,  6,  8,  and  i|t. 

See  Figure  2k  (Reference  1)  for  the  directions  of  the  beam  axes. 

3.  The  ub- at  the  time  of  rupture.  This  is  the  total  deflection, 

XJX 

in  beam  axes,  in  the  direction  of  the  rupture. 
k.  Uie  Umax- •  which  defines  the  maximum  allowable  beam  deflection 

jo1 

in  the  l1,11  direction.  This  is  an  input  constant. 


If  a  control  volume  mass  penetration  occurs,  the  mass  which  penetrates  and 
the  time  are  printed  out.  At  the  end  of  the  run,  the  ruptures  and  mass 
penetrations  that  occurred  during  the  run  are  summarized  in  tables  for 
ready  reference. 

Sample  Cases 

Two  sample  cases  are  presented.  The  first  sample  case  represents  a  com¬ 
bined  vertical  and  lateral  impact  condition.  The  vertical  impact  velocity 
is  23  ft/sec  and  the  lateral  impact  velocity  is  l8.6  ft/sec.  The  mathe¬ 
matical  model  is  the  31  mass  representation  shown  in  Volume  I,  Figure  6 
and  Table  VT.  The  sample  output  for  one  time  period  is  shown  in  Figure  85. 
The  input  data,  which  follows  the  User's  Guide  Input  format  (Figure  84),  is 
presented  in  Figure  86. 

Sample  output  plots  are  provided  for  the  engine  mass  vertical  velocity, 
engine  mass  vertical  acceleration  and  engine  mount  support  vertical  deflec¬ 
tion  in  Figures  87,  88  and  89,  respectively.  Figures  90)  91  and  $2  provide 
similar  sample  plots  for  the  lateral  velocity,  lateral  acceleration  and 
mount  lateral  displacement,  respectively.  Figures  93  and  94  show  the  for¬ 
ward  and  aft  DRI  plots,  respectively,  for  the  comoined  vertical  and  lateral 
impact  velocity  case.  The  engine  vertical  responses  shown  in  Figures  87 
(velocity),  88  (acceleration)  and  89  (displacement)  can  be  compared  to 
Figures  17j  18  and  19  in  Volume  II,  Reference  1  to  show  that  the  revised 
program  KRASH  performs  in  the  same  manner  as  the  final  submitted  program 
KRASH  from  the  previous  study  (Reference  l) . 

The  second  sample  case  provided  in  this  section  is  the  three-dimensional 
velocity  impact  condition.  The  initial  c.g.  velocity  for  this  case  is  a 
combined  40  ft/sec  velocity  in  which  the  velocity  vector  components  are: 

Vx  =  27  ft /sec  V  =  l8.6  ft /sec  V^  =  23  ft /sec 

The  initial  vehicle  attitude  is  flat  (zero  roll,  pitch  and  yaw  angles),  and 
the  initial  angular  velocities  are  zero.  Figure  101  (page  1 U 0)  in  Volume  I 
shows  the  mathematical  model  representation  for  this  condition  The  model 
consists  of  32  lumped  masses.  It  is  the  same  model  used  in  the  previous 
sample  problem  except  for  the  extra  mass  location  which  represents  the 
rotor  blade  during  a  potential  mass  penetration  (of  the  upper  cabin)  con¬ 
dition.  The  input  data  for  this  sample  problem  is  provided  in  Figure  95* 

A  sample  output  for  one  time  period  is  shown  in  Figure  96.  Figures  97 
through  99  show  the  engine  mass  vertical,  lateral  and  longitudinal  veloc¬ 
ities,  respectively.  Figures  100  through  102  present  the  engine  mass  ver¬ 
tical,  lateral  and  longitudinal  accelerations,  respectively.  Figures  103 
through  105  show  the  engine  support  mount  vertical,  lateral  and  longitudinal 
deflections,  respectively.  The  forward  and  aft  DRI's  are  presented  in  Fig¬ 
ures  106  and  107. 

Program  Listing 

A  complete  listing  of  the  computer  program  KRASH,  including  subroutines, 

'is  presented  in  Figure  108, 
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Figure  85.  (Continued) 
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Figure  85.  (Continued) 
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Figure  85.  (Continued') 
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Figure  85.  (Continued) 


Combined  Vertical  and  Lateral  Impact  Sample  Case . 
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Figure  86 .  ( Continued ) . 
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Figure  <?k.  Af't  DRI  Time  History  Plot,  Combined  Vertical 
and  Lateral  Impact  Sample  Case. 
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Figure  95*  Input  Data,  Combined  Vertical,  Lateral  and 
Longitudinal  Impact  Sample  Case. 
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Figure  95.  (Continued) 


f  I 


OOH  O  *| 


090  OOi 
•  •  •  •  4 

O  •  O  OOi 

f 

o 
o 

1 

*oo  o  o) 
(XOO  O  oj 


o«o 
•  *  • 
0*0 


O  Of 
•  •  • 
oom 
I 


o+\ 


o  c  o  •  o 
•  ••#••• 
o  o  O  r*  *  O 


Ooo  • 
•  *  •  • 
ooo« 


o' 

8; 


••ooo 

•  t  «  I 
NOflO 


o«oo 


8  & 


•  •  •  • 

WON6  OjOOOOO 


OOO  9 

•  •  •  i 

ooo« 

1 1 


oooo  « 

•  lilt 

090  0  4 


O  ooo  o 

•  •  I  •  « 

ooooc 


09600 

!••••• 

ooooo 


o  o«o  o 
•  •  •  •  • 
OONOO 


ooooo 
•  *  •  »  • 
ooooo 


04000$ 
.«*••••  4 
ONOOOH 


o  e  4  o  a  d 
0040-d 


booMOd 

•  •  •  •  •  4 
OOO  NOoj 

s  s 

§  i 

,  3  i 

lo  o  •  o  *  a 

•  ••••• 

oo-o^o. 


ooooo* 

OhOOON 


oooooo 

•  •••••I 

••ooooo 


onood 

•  •  •  •  4 

OlAOOOj 


t 

i  t 


is  l 

009  09 
•  •  •  •  • 
0040  f* 


;s  a 

'o  d 

§  i 

o  m\ 
0  0  90  91 
•  *  •  •  *i 
OOMO  w 


0*0  0  o\*i 

0*00  0,(0 


40090 
«*••••• 
«f»  0.0  c  o 


oto, 


ooo  « 

•  •  •  4 

oo  *  c] 


ooo  S 
•  •  •  4 
oood 


i  2 


0*4  o«o{ 

•  •  •  •  4 

on  o|q 


i 

8 


H  o  OOd 
•  •  ♦  •  4 

NO  '©0*t 


ooo] 


1  t 

_  _  _  j  «  o  oo* 
•  ••• 
CONO-OO  OON 


I 


8 

o 

f 


0*0  0:0  « 
•  •••!•• 
ONO  0,0  * 
I 


o 

8 

o 

N  o « 

-  4  •  •  , 

o  *-•  *  o  o  o 


0*0 
•  *  ♦ 
ONO 


'8 

o 


>  £  o  o 

00*4400 


004 
•  •  • 
OON 


O  MO 

f  •  • 

ONO 


9  OO 
•  •  • 
mop 


00  9  0  4 

009  0*4 


ooo 
•  ••- 
O  O  ON*  0\ 


oooo 

•  I  •  • 

0  0  0*1 


oooo 

•  •  •  • 


oooo 
•  •  *  • 
000)0 


90  9  0 

•  «  •  I 

•40  NO 


oooo 
•  •  •  • 
00941 
4 


oooo 

*  •  * 

O  OOONO) 


OOOOO 

4IOOQH 


9 

ooooo 
•  •  »  •  • 
O  N  ©  «  o 


ooooo 

i  •  •  *  •  • 
OOOOO 


s  It 

o  jo 

8  8 


NO  OOO 

I  9 


ooooo  , 

0  0  0  0  0-4 


o  m  do  o  *J 

•  •  •  •  •  4 
OOOO  ON 


0004)00' 
•  •••••, 
OO  Of  oo 


s  *  , 

oo  f or  o 
•  ••••• 
OOOOHO 


J 


pmooom 

t  I  •  9  •  I 

O  **  OOO  m 


f  o  ooo  o 
•  ••»»• 
>40  0,00  0. 


143 


Figure  95.  (Continued) 
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Figure  96.  (Continued). 
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Figure  96.  (Continued) 
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Figure  97.  Engine  Mass  Vertical  Velocity,  Combined  Vertical 
Lateral  and  Longitudinal  Impact  Sample  Case . 
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Figure  103 .  Engine  Mount  Vertical  reflection.  Combined  Vertical 
Lateral  and  Longitudinal  Impact  Sample  Case. 
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Figure  lOU.  Engine  Mount  Lateral  Deflection,  Combined  Vertical 
Lateral  and  Longitudinal  Impact  Sample  Case. 
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Figure  107.  Aft  DRI,  Combined  Vertical,  Lateral  and 
Longitudinal  Impact  Sample  Case. 
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Figure  108.  (Continued) 
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Figure  108.  (Continued) 
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Figure  108.  (Continued) 
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Figure  108.  (Continued) 


Figure  108.  (Continued) 
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area;  ’for  stiffened  panel,  area  of  stiffener  plus  area 
of  sheet  corresponding  to  stiffener  spacing,  in. 

plate  length,  in. 

2 

area  of  stiffener,  in. 

2 

area  of  web,  in. 

width  of  skin,  in. 

width  of  stiffener  web,  in. 

constant  (=  3*0) 
number  of  cuts 
rivet  diameter,  in. 
modulus  of  elasticity,  psi 
distance  to  neutral  axis,  in. 

number  of  flanges 
crush  strength,  lb 

effective  rivet  off-set,  in. 

acceleration  level 

2 

gravity  term  =  386.2  in. /sec 
distance,  in. 

4 

moment  of  inertia,  in. 

4 

moment  of  inertia  about  c.g.,  in. 

eccentricity  factor 
kinetic  energy,  in. -lb 
wrinkling  coefficient 

one -half  column  length,  in. 
effective  length  =  L/  n/c",  in. 
effective  moment,  in. -lb 
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plastic  hinge  moment,  in. -lb 
neutral  axis 

effective  axial  force,  lb 
plastic  hinge  axial  force,  lb. 

rivet  pitch,  in. 
applied  axial  force,  lb 
monolithic  failure  load,  lb 

wrinkling  failure  load,  lb 
distributed  load,  Ib/in. 

stopping  distance,  core  radius  of  given  cross  section,  in. 

distance  from  midcross-section  of  skin  to  neutral  axis 
of  stiffener,  in. 

ratio  of  a/ a 
'  cy 

ratio  of  critical  stress  (a  )  to  yield  stress  (o  ) 

cr  cy 

minimum  core  thickness,  in. 

skin  thickness,  in. 

stiffener  thickness,  in. 

impact  velocity,  in. /sec 
weight,  lb 

effective  half  width  of  skin,  in. 
lateral  deflection 

function  of  column  length  and  end  shortening,  in. 
defined  in  equation  51 
plasticity  reduction  factor 
cladding  reduction  factor 

Pi  =  3*lU 
tan”1'  Z 


radius  of  gyration  =  in* 


partial  derivative 


SUBSTRUCTURE  ANALYSIS  (Continued) 


6u  incremental  displacement 

6,  A,  u  displacement,  end  shortening 

v  Poisson's  ratio 

x  distance  to  c.g.,  in. 

o  compressive  stress,  psi 

o k  stress  at  outermost  fiber  of  sheet  skin,  psi 

o  column  buckling  stress,  psi 

CO 

0cr  buckling  stress,  psi 


compressive  yield  stress,  psi 
effective  compressive  yield  stress,  psi 
effective  compressive  yield  stress  for  stiffener  web,  psi 
effective  compressive  yield  stress  for  skin,  psi 
Euler  stress,  psi 

crippling  (or  failure)  stress,  psi 
failure  stress  of  riveted  panel,  psi 
failure  stress  due  to  wrinkling  mode,  psi 
buckling  stress  at  L'/p  -  20 


Subscripts: 


b 

f 

i 


bending 

filler 

.th  . 

l  segment 


max. 

min. 

P 

s 

w 


maximum 
minimum 
plastic  state 
skin  or  sheet 
stiffener  web 


M 


SUBSTRUCTURE  ANALYSIS  (Continued) 

Constants  (function  of  stiffened  panel  configuration) 


"Sci’  Vi’  *zi 

III 
xyi’  yzi*  zxi 


rotational  transformation  matrix  from  body  axes  to 
ground  axes 

transpose  of  matrix 

total  crash  spring  (external  spring)  energy  absorbed 
total  damping  energy  dissipated 
dynamic  response  index 
total  system  energy 

til 

internal  beam  damping  force  (or  moment)  for  the  ij 
beam  in  the  0  th  direction 

crash  spring  forces;  spring  ij  in  the  k  direction 

force  (or  moment)  at  point  j  due  to  beam  ij,  in  the 
0th  direction 

"til 

moments  of  inertia  of  lumped  mass  m.  about  1  1  body 
fixed  axes  x 

«j. 

product  of  inertia  of  lumped  mass  m.  about  i'  body 
axes 

total  kinetic  energy 

length  of  vector  from  m.  to  ground  contact  point  C' 

,th  ,  , 

l  lumped  mass 

number  of  lumped  masses 

total  potential  energy 

total  strain  energy  absorbed 

term  in  expression  for  crash  spring  energy 
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W. 


,th 


weight  of  i  lumped  mass 


x 


v  y,,  z 


U^,  V,,  W 


1' 

rv 


V  V  ri 


Ainp^,  Ainq^  Ain^ 

Avi«i 


ground  coordinates  of  rru 
til 

i  body  axes  component  of  absolute  translational 
velocity  vector  of  mass  i 
til 

i  body  axes  components  of  absolute  angular 
velocity  vector  of  mass  i 

til 

crash  (external)  forces  and  moments,  i  body  axes 
=  Iavc. 

j  i|  th 

incremental  displacement  of  point  i,  i  n  body  axes 

incremental  rotation  of  point  i,  i  body  axes 

incremental  displacement  vector  of  point  j  with 
respect  to  point  i,  due  to  deformation  of  beam  ij, 
in  the  0th  direction 


$ 
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Bj_,  Ci 

[Ail 

KJ 

W 

t*I] 

[»•] 

[3'] 


ik 

i 

lk 

dvciJk 

d™ijk 

M 

b’3 

^ijkl 

^ijkl 

PSPijk 

rePOik 


Terns  used  in  Euler's  equations  of  motion 

Rotation  transformation  matrix  from  body  axes  to 
ground  axes 

Matrix  relating  (0A,  to  (p1#  qi#  r,)  in 

equation  (92)  1 

Time  derivative  of  M 

Rotation  transformation  matrix  from  ith  body  axes 

to  c.g.  axes 

Rotation  transformation  matrix  from  beam  ij  axes  to 
body  i  axes 

Rotation  transformation  matrix  from  c.g.  axes  to 
ground  axes 

Matrix  relating  (0*,  ■*',  *//')  to  (p',  q',  r')  in 
equation  (88) 

End  point  of  kth  spring  on  ith  mass 

Ground  contact  point  of  kth  spring  on  ith  mass 

Ground  axes  components  of  rector  from  mi  to 

Ground  axes  components  of  vector  for  to  Clk 

Derivative  matrix 

Derivative  matrix 

Running  time  sum  of  AfMijj^ 

Value  of  FMjjjq  at  time  of  loading  reversal 

Body  i  axes  components  of  spring  force  at  ground 
contact  point 


Axial  compressive  force  in  kth  spring  on  ith  mass 
Value  of  FSPOjv  at  time  of  loading  reversal 


USER'S  GUIDE  AND  REFERENCE  1  TERMS 


1 


reKVlK 

G 

H 

Hexi>  Hezi 

Xxi>  *yi>  Tzi 

Ixyi»  Xzxi 

keik 

fcj 

CkriJ 

1ik 

•^ci*  ^yi*  ^i^ik) 

^ci 

LIFTi 

“l 

HlUik 

N 

%k 

nxj  ^y)  nz 

0 


Final  value  of  FSPO^k  in  input  table  of  Sj_k  vs.  FSPO^k 

Center-of-gravity  of  total  vehicle 

Origin  of  helicopter  coordinate  system  (F.S.O, 

B.L.O,  W.L.O) 

Angular  momenta  of  mi  due  to  rotation  of  masses  in¬ 
ternal  to  mj 

Moments  of  inertia  of  lumped  mass  mi,  about  ith 
body  fixed  axes 

Products  of  inertia  of  lumped  mass  m*,  about  ith 
body  fixed  axes 

Linear  unloading  stiffness  for  kth  spring 

Six  by  six  linear  stiffness  .matrix  for  beam  ij 

Six  by  six  diagonal  stiffness  reduction  matrix  for 
beam  ij 

I 

Length  of  vector  from  to  ground  contact  point  Cik 
Free  length  of  kth  spring  on  ith  mass 
Aerodynamic  lift  constant 

Aerodynamic  lift  on  14,  positive  up,  in  ground  axes 
ith  lumped  mass 

Ground-spring  friction  coefficient  for  kth  spring 
on  ith  mass 

Total  number  of  lumped  masses 

Unit  vector  triad  fixed  in  ith  body  coordinate  system 
Unit  vector  triad  fixed  in  ground  coordinate  system 
Origin  of  ground  coordinate  system 
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Pi#  *i#  ri 


ith  body  axes  components  of  absolute  angular  velo- 
city  vector  of  mass  i 


P*#  d'#  r* 

•a 


c.g.  axes  components  of  initial(t*0)  vehicle  angular 
velocity  vector 

Contact  point  velocity  matrix  used  in  equation  (60) 

Axial  external  spring  compression,  kth  spring  on 
ith  mass 


8ik 


Value  of  sik  at  time  of  loading  reversal 


•ik 


Final  value  of  in  input  table  of  ®lk  vs.  FSPOjj^ 

kth  spring  axial  compression  measured  relative  to 
current  load  stroke  curve  origin 


Horizontal  shift  of  coordinates  with  respect  to 
Sik  coordinates 


t 


Time 


[TuJ 

«i#  vi#  wi 


Static  balance  matrix  used  in  equation  (30b) 

Body  1  axes  components  of  absolute  translational 
velocity  vector  of  point 


*>1J 

*>ijl 

~ik 

vtijl 

^ijl 

vcijk 


*i#  y±,  zj 

Running  time  sum  of  Avbjj 

Value  of  vbjji  at  time  of  loading  reversal 

Magnitude  of  ground  plane  contact  point  velocity 

1th  total  beam  deflection  measured  relative  to 
current  load- stroke  curve  origin 

Horizontal  shift  of  vh^  coordinates  vith  re¬ 
spect  to  vb^jjL  coordinates 

Ground  coordinates  of  point 


9 
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vo»uk 

*YC' 

oyc> 

Of* 


ybijf 

*G#  yo#  ZG  (v«j) 

*G>  Yq,  ZG 

i»  it  ii  .  » 

*G>  Yq$  *G  (v«PPij) 

xi#  yi>  zi  (v»ij) 
*i>  Ylt  H  (viPij) 

*i#  y£>  zi 

«  «  i 

xij*  Vit*  zij 

/x{j*  ziA 

\Li i*  Mij»  hj 

(Ai>  Z'A 

\Lji>  Mji,  V'J 


Ground  axes  components  f  of  absolute  velocity  of 
ground  contact  point 

Velocity  vector  of  cifc  vith  respect  to  »i 
Velocity  vector  of  cjk  vith  respect  to  ground 
Velocity  vector  of  mj,  vith  respect  to  ground 
Weight  of  1th  lumped  mass 
Total  vehicle  veight 
Beam  ij  coordinates 

Ground  coordinates  of  initial  (t  »  o)  c.g.  position 

Ground  axes  components  of  initial  (t  ■  o)  c.g. 
velocity  vector 

Helicopter  axes  coordinates  of  vehicle  c.g. 

(point  G) 

Ground  coordinates  of  a* 

Coordinates  of  m*  in  center-of- gravity  coordinate 
system 

Coordinates  of  in  helicopter  coordinate  systen 

Ground  coordinates  of  vector  from  point  i  to 
point  J 

ith  body  coordinates  of  vector  from  point  i  to 
point  J 

Total  (summed  over  time)  internal  forces  and  moments 
at  point  i  due  to  beam  i  j ,  i^h  body  axes 

Total  (summed  over  time)  internal  forces  and  moments 
at  point  j  due  to  beam  i j ,  jth  body  axes 
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/Xi>  Yi>  Zi\ 
\Li>  N i  / 

XAi  YAi>  2Ai 

/xci*  Yc u  zci\ 
V  LCi>  MCi>  NCi/ 


XGl>  V  ZGi 


(*n>  Yii'  ziA 

\Ln>  Nxi/ 


Total  forces  and  moments  on  mass  i,  in  ith  body  axes 


Aerodynamic  forces,  i  body  axes 

Crash  (external)  forces  and  moments,  i  body  axes 


Gravity  forces,  i  body  axes 

Internal  forces  and  moments,  i  body  axes 

Ground  axes  components  of  spring  force  at  ground 
contact  point  c’k,  positive  up,  left  and  aft 

Vertical  distance  from  c.g.  to  lowest 

Determinate  expression  used  in  equation  (68) 

Incremental  forces  and  moments  at  point  j  due  to 
beam  ij 

kth  incremental  load  due  to  1th  incremental  deflec¬ 
tion  for  beam  ij 


A^b^,  A^b^,  j  Incremental  rotations  of  point  j  with  respect  to 

point  i,  in  beam  ij  axes 

Mi  Incremental  change  in  ith  mass  Euler  angles 

At  Numerical  integration  time  interval 


^io 


Incremental  change  in  ith  mass  Euler  angles 

Numerical  integration  time  interval 

Six  element  vector  made  up  of  Axb.  ,,  Ayb.  ,,Azb. 

Aifey,  d«b1J(  1J  1J  lj’ 

Incremental  displacement  vector  of  point  j  with 
respect  to  point  i,  due  to  deformation  of  beam  ij 

Incremental  displacement  vector  of  point  j  with 
respect  to  point  i 
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A^ij 

Axb^Aybij^zb^ 

Incremental  displacement  vector  of  point  J  with 
respect  to  point  i,  due  to  rotation  of  mass  i 

Coordinates  of  Avbitj  in  "beam  ij  axes 

Axi#  4yi#  Az± 

Incremental  displacement  of  point  i,  ground  axes 

Axu'Ayij'Ai!ij 

Incremental  displacement  of  point  j  with  respect 
to  point  i  in  ground  axes 

l  *h  j>  ^ij 

Incremental  internal  forces  and  moments  at  point 

J  due  to  beam  ij,  in  beam  ij  axes  (elements  of 
vector) 

ULij<AKiJ'AN?j/ 

Incremental  internal  forces  and  moments  at  point  j 
due  to  beam  ij,  ground  axes 

/SSfy^TS^JN 

AM?j  A"?/ 

Incremental  internal  forces  and  moments  at  point  i 
due  to  beam  ij  ,grcund  axes 

/OxJj.AyJj.AzJA 
(  Ai4j,AMij,ANi/ 

Incremental  internal  forces  and  moments  at  point  i 
due  to  beam  i j ,  ith  body  axes 

/ AXji,  AYjj,  AZjA 
ANjJ 

Incremental  internal  forces  and  moments  at  point  j 
due  to  beam  i j ,  jth  body  axes 

01,  •%,')& 

Euler  angles  from  ground  axes  to  body  axes  (time 
varying) 

^tJ*  "*« 'fij 

Euler  angles  from  ith  body  axes  to  beam  ij  axes 
(constant) 

t',  ’ 

Euler  angles  from  ground  axes  to  c.g.  axes  (con¬ 
stant);  initial  (t=o  )•  attitude  of  vehicle 

&  fi 

Euler  angles  from  c.g.  axes  to  ith  body  axes  (con¬ 
stant  ) 
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